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MULTIWIRE PROPORTIONAL CHAMBERS FOR TRAJECTORY 
MEASUREMENTS IN HIGH ENERGY COSMIC RAY EXPERIMENTS 


1 . INTRODUCTION 

One of the most important measurements in high energy physics experiments 
- be it in elementary particle physics or in cosmic ray physics - is the 
exact determination of the trajectory of the particle or particles under 
investigation. In high energy cosmic ray experiments the trajectory of 
particles has to be known for the following reasons: 

1. In order to determine the charge of a high energy (i.e., relativistic) 
charged particle with high resolution, it is necessary to know the 
location and length of the path of that particle through the charge 
measuring detectors, which can be scintillators, Cerenkov counters 

or ionization chambers. 

2. In a magnetic spectrometer experiment such as the Superconducting 
Magnetic Spectrometer Experiment BCR-S^on HEAO-B the exact incoming 
and outgoing trajectory of the particle has to be known in order to 
determine the "bending" of the particle's path in the magnetic field. 
The amount of bending is a measure of the rigidity of the traversing 
particle and the direction of bending allows to distinguish between 
positively and negatively charged particles. 

2 

3. In an ionization spectrometer experiment such as the ACR-6 experiment 
on HEAO-A the trajectory of the incoming particle has to be measured 
in order to determine the shower axis of the resulting nuclear-electro 
magnetic cascade in the ionization spectrometer. The knowledge of the 
shower axis allows for certain corrections of the energy measurement 
due to the leakage of the cascade out of the side and bottom of the 
spectrometer. 

4. In order to recognize and to determine the nature of simultaneous back 
ground particles (e.g., from backscattering from heavy absorbers like 
the Total Absorbtion Shower Counter in BCR-5 or the Ionization Spectro 
meter in ACR-6), it is necessary to simultaneously measure the tra- 
jectory of more than one particle which will then help to determine 
the trajectory of the primary particle. 
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A very important reason for measuring the trajectory of high energy 
cosmic ray particles is the investigation of anisotropy of the cosmic 
ray flux. Without good directional resolution it will not be possible 
to investigate anisotropies of a scale size that is smaller than the 
angular aperture of the experiment. 

6. In magnetic spectrometer experiments like the one that was developed 
3 

by R. Golden of NASA's Manned Space Center in Houston for balloon 
flights and which will be developed for shuttle flights, spatial de- 
tectors are used in combination with nuclear emulsion to greatly in- 
crease the resolving power of such spectrometers and to extend their 
range to higher energies. The spatial detector locates the particle 
with sufficient accuracy to enable automatic location of tracks in 
the nuclear emulsion which has superior spatial resolution. 

In the past mostly spark chambers were used to measure the trajectory of 
high energy particles and it is interesting to see that the original HEAO 
proposals for the BCR-5 and ACR-6 experiments contained spark chamber 
systems for trajectory measurements. 

Since the successful operation of large area multiwire proportional 
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chambers (MWPC) by Charpak et al. at CERN in 1968, these detectors have 
come into widespread use, and a significant amount of development work has 
been done in recent years. The work at TRW was directed towards the de- 
velopment of large area MWPC-systems with high spatial resolution for 
eventual use in balloon and satellite experiments like the ones that were 
proposed for the HEAO missions and which will later be flown on shuttle. 
Much of the work on high resolution MWPC's is also applicable to the de- 
velopment of imaging proportional chambers for x-ray experiments in x-ray 
astronomy. Especially the delay-line readout technique for MWPC's is 
ideal for x-ray imaging proportional chambers. 

This report describes the work that was performed by TRW under NASA con- 
tract No. NAS9-11315, towards the development of high resolution, large 
area MWPC's. In the. following sections we present first a set of typical 
requirements for a spatial detector system for a cosmic ray experiment 
(Section 2). We discuss the advantages of MWPC's over spark chambers as 
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spatial detectors. 

In Section 3 we give a short description of the principle of operation 
of MWPC's, and we present a summary of various readout systems which 
have been used and/or have been proposed. We close with a brief char- 
acterization of the lumped element delay-line concept which has been 
used at TRW. This leads into a summary of the parameters which influence 
the achievable spatial resolution. 

After that we describe in more detail the experimental and theoretical 
work that has been carried out at TRW. Gas gain calculations and satur- 
ation effects due to space charge buildup are presented in Section 4, and 
the results of our gas gain measurements are presented. We then discuss 
in more detail limits on spatial resolution and the results of spatial 
resolution measurements (Section 5 and 7) . A complete analysis of the 
lumped element dealy-line is given in Section 6. Section 8 presents our 
attempts to evaluate the delta-ray problem with the help of Monte Carlo 
calculations which were used to simulate delta-ray effects. In Section 
9 we describe the results of a set of measurements we have carried out 
with MWPC's in a high energy (2.1 GeV/Nucleon) , heavy ion beam (C 12 and 
016) at the Bevatron. Some design considerations and results of load 
calculations for MWPC frames are presented in Section 10. All references 
which were used during preparation of this report are also listed. In 
addition we have compiled a bibliography of selected papers dealing with 
various aspects of MWPC's. 

2. REQUIREMENTS ON SPATIAL DETECTORS IN SATELLITE EXPERIMENTS 

The use of spatial detectors in high energy cosmic ray experiments on 
satellites imposes strict requirements on the design, construction and 
performance of the detectors due to limitations in available weight, vol- 
ume and power and due to the requirement of high reliability and long 
lifetime. 

A representative set of requirements are the requirements for the spatial 
detectors in the Superconducting Magnetic Spectrometer Experiment (BCR-5) 
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which was selected to fly on the HEAO-B spacecraft. 

The experiment required 8 sets of detectors to measure the position of 
cosmic rays with unit charge up to Z=26 with a resolution of 0.1 mm in 
each set in X-and Y-direction. Assuming the position can be determined 
with an accuracy of 0.2 mm in a single spatial detector gap each set has 
to contain 8 gaps (4-x and 4-y) to achieve the required resolution per 
spatial detector set. This leads to a total of 64 gaps with the fol- 
lowing parameters: 


Effective area: 

Maximum Thickness: 
Weight: 

Average weight per gap: 
Power dissipation: 


4 sets 66 cm x 53.3 cm 
4 sets 45 cm x 42 cm 
4 inch per set 
52 lbs for large sets 
48 lbs for small sets 
approx. 6 lbs 
4 watts per set 
or 500 mW per gap 


It was further required that the septum width of the detector frame along 
the cryostat of the superconducting magnet be kept as small as possible 
in order to minimize the loss of high magnetic field regions for pene- 
trating particles. The electronic readout system was required to measure 
two simultaneous trajectories for Z 2 events. The lifetime of the ex- 
periment in orbit was to be at least one year and preferably two years. 


2 . 1 MWPC vs Spark Chambers in Cosmic Ray Experiments 

In this chapter we discuss the main advantages of MWPC's over spark 
chambers in high energy cosmic ray experiments (balloon-borne or satellite- 
borne) . Because of these advantages the decision was made to replace spark 
chambers with MWPC's in the BCR-5 and ACR-6 experiments on HEAO. In ad- 
dition to the advantages discussed here there are others which are re- 
sponsible for the usefulness of MWPC's in nuclear and elementary particle 
physics experiments at accelerators, like the considerably higher event 
rate that a MWPC can process. 

MWPC's operate with a dc-potential between anode and cathode and are con- 
stantly sensitive to penetrating particles. Because the recovery time is 
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short (about 100 nsec), dc-operation is possible even in the presense 
of high background rates, and other detectors within the experiment can 
decide whether or not the trajectory of a penetrating particle is to be 
read out . 

Spark chambers, on the other hand, are triggered devices and in order to 
sensitize the spark chamber, a high voltage pulse has to be applied im- 
mediately after the passage of a particle whose trajectory has to be de- 
termined. This not only requires fast, high voltage trigger circuits with 
very high reliability (now under development), but it also creates a very 
severe electromagnetic interference (EMI) problem. EMI from the spark 
chambers requires special precautions to be taken to protect other elec- 
tronic circuitry and to reduce errors due to EMI in the measurement of 
low-level analog signals. Spark chambers have to be either extremely 
carefully shielded or the measurement and analysis of low-level analog 
signals from detectors like scintillation-counters has to be completed 
before the spark chambers are triggered. This requires fast electronic 
circuitry. While shielding adds weight and complexity, fast electronics 
consumes considerable amounts of electrical power. Both might be pro- 
hibitive in a given experiment. 

MWPC's only require stable and programmable high voltage power supplies, 
which are readily available, and they do not emit EMI which allows the 
use of low power, slow electronics for the other detectors in the same 
experiment which leads to lower power consumption. The complexity of the 
readout electronics is roughly the same for both types of detectors. 

Another advantage of MWPC's is their expected longer lifetime. Spark 
chambers use small amounts of organic quenching agents like ethanol for 
stable operation. These additives, however } can be disassociated during 
spark breakdown and have to be replaced. The disassociation products also 
can be the source of serious contamination through the buildup of perman- 
ent deposits on spark chamber wires which leads to unstable operation. 
MWPC's can successfully be operated with gas mixtures of argon, xenon, and 
carbon dioxide, i.e., without organic quenchers. These gas mixtures are 
not very prone to decomposition and it is believed that even sealed, large 
area MWPC's can be constructed and operated over long periods of time 
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without significant deterioration. It should be pointed out, however, 
that decomposition and contamination of MWPC gases can present a serious 
problem if gas mixtures with complex organic components like methane, 
isobutane, or freon are used, especially if the chambers have to be oper- 
ated in a high gas gain mode. 

A further very important advantage of MWPC's is their better ability to 
discriminate against delta-rays. Delta-rays are energetic electrons 
which are created when charged particles traverse any kind of material. 

The number of delta-rays increases with (Z is the charge of the pene- 
trating particle) which means that especially high-Z cosmic ray nuclei can 
be accompanied by considerable numbers of them. The task of the spatial 
detectors is then to determine the position information of the high-Z 
particle out of the background of simultaneous delta-rays. Normally 
the difference in energy loss between a delta-ray and a high-Z particle 
is quite significant in the sensitive volume around a MWPC anode wire, 
and because of the proportionality between energy loss and MWPC output 
signal it is in principle possible to discriminate against delta-rays by 
electronic means and to measure the trajectory of a high-Z cosmic ray 
nuclei unambigiously. Delta-ray discrimination with spark chambers is 
inherently much more difficult because the spark mode of operation produces 
an output signal of a spark chamber that is not proportional to the energy 
loss in the spark chamber gap. 

3. MWPC - PRINCIPLE OF OPERATION 

Single wire proportional counters have been used for a long time and 
their operation is well understood. The gas multiplication process around 
the anode wires of a MWPC is the same as in a single wire cylindrical 
counter and will, therefore, not be discussed in any detail, unless it is 
necessary for the understanding of a MWPC as a spatial detector. Charpak^ 
has written a very good review article about MWPC's and their application 
in high energy physics experiments, and we will frequently refer to this 
article . 

The typical MWPC configuration that has been under investigation at TRW is 
shown in Figure 3-1. The anode and cathode wires were wound on round fiber- 
glass frames and soldered to a pc-board, which was glued into the fiber- 
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ANODE 
WIRE PLANE 



ANODE 
WIRE PLANE 



STANDARD LABORATORY CHAMBER: 


ACT IVE AREA 
L 
S 


2" X 2" 
0.5 CM 
0.4 CM 


d : 20 MICRONS 


FIGURE 3-1 MWPC- Con figuration 
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glass frames. The fiberglass rings were then clamped together with fi- 
rings in between to seal off the chamber volume. Two additional fiber- 
glass rings with thin Mylar films were used as gas windows. The chambers 
were operated in the gas flow-through mode. Typical chamber parameters 
are : 

Anode - Cathode Gap: 5 mm 

Anode wire spacing: 4 mm 

Cathode wire spacing: 2 mm 
Anode wires perpendicular to Cathode wires 
Anode wires: 20y diameter gold plated tungsten, 

the two end wires on each side of 
the anode plane wire 50y diameter 
wires 

Cathode wires: 50y diameter gold plated tungsten 

Also shown is the field configuration around the anode wires. Although 
the anode wires are not separated by cathodes it can be seem that the 
electric field forms a well defined sensitive cell around each anode wire. 
It is also indicated that the electric field close to the anode wire is 
nearly cylindrical. 

When a charged particle enters the MWPC gap it produces a number of 
electron-ion pairs through ionization of the chamber gas. The electrons 
then start drifting towards the anode and the positive ions towards the 
cathodes. Because of the high field strength close to the anode wire the 
primary electrons are accelerated to sufficient energies to ionize further 
gas atoms and to initiate a cascade process which greatly multiplies the 
original number of ion-electron pairs and, therefore, leads to an observ- 
able signal. The cascade process occurs in very close proximity to the 
anode wire, i.e. in a cylinder around the wire with a radius of about one 
or two wire diameters. (See Ref. 5). Direct collection at the anode of 
the electrons is only responsible for a small fraction of the total pro- 
portional chamber signal, the largest part of the signal is due to the 
movement of the positive ions along electric field lines away from the 
anode. This can be explained as follows: The movement of electrons and 

positive ions induces currents (1) in the anode wire which collects the 
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electrons, (2) in neighboring anode wires in a MWPC, and (3) in both 
cathode planes. The induced current is 

• - - AQ. Ar Ne AIJ 

1 - At ~ Ar At - V 0 Ar W 

Ne is the charge moving through a potential drop AU around the anode wire 
inducing the charge AQ which varies with time and causes the current i to 
flow. V Q is the applied voltage and w is the mobility of the charges. 

In low-field regions the mobility of the positive ions w + is much lower 
than the mobility of the electrons w_. However, in the high field region 
close to the collecting anode wire the mobilities are about equal: w + « 
w_. The electrons of the multiplication region see only a small value of 
AU before being collected while the positive ions must move through the 
full potential drop across the chamber. The drift of the positive ions 
away from the anode is, therefore, responsible for most of the total MWPC 
signal. Since the electrons and initially also the positive ions move 
very fast one sees a fast rising current pulse with a slow decay. The 
slow decay represents the movement of the positive ion cloud in low field 
regions. Figure 3-2 shows some current pulses taken from one of our MWPC's. 
The pulses were produced by 5.9 KeV x-rays from a Fe-55 source. The 
chamber gas was 90% Ar + 10% CO 2 and operating voltage was 2100 volts. 

Table III-l shows the time development of the charge pulse in our chamber 
configuration, based on equations given by Charpak 5 . The important thing 
to notice here is the fact that it takes many usee to collect the total 
charge and that only about 40% of the charge is collected during the first 
usee, i.e. during the fast part of the pulse buildup. 

In many MWPC applications it is not feasible to make use of the total 
charge. Typically, the output of a charge sensitive preamplifier is 
shaped by a post amplifier with RC shaping time constants of one or two 
microseconds in order to operate at reasonable pulse rates. This means 
that normally less than 10% - 50% of the charge is seen, depending on 
chambers configuration and operating parameters. This is also important 
with respect to the term "gas gain" of a proportional counter. Gas gain 
is the numerical factor which states by how many times the number of pri- 
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TABLE III-l 

Time development of the MWPC charge pulse Q(t) in percent of the total 
charge (see also Ref. 5) 
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> 

o 
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FIGURE 3-2. Photograph of Fe-55 current pulses in 
90% Ar + 10% C0 2 at 2100 Volts 
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many electrons has been multiplied in the cascade process. Numbers which 
are quoted often do not represent the total gas gain, but rather a rela- 
tive gas gain referring. only to the fast part of the charge pulse because 
shaping amplifiers are used. 

The movement of the positive ions away from the collecting anode wire in- 
duces a negative current pulse in that wire and positive current pulses 
in the neighboring anode wires and the two cathode planes (total induced 
charge = 0) . The way the positive charge is divided between the cathode 
planes and the neighboring anodes depends on the geometrical configur- 
ation of the MWPC. In order to have most of the charge induced in the 
cathode planes it is important to make the ratio of anode wire spacing 
to anode-cathode gap as large as feasible. 

The output pulse of a MWPC is over a large region strictly proportional 
to the initial ionization. However, if the charge density of the positive 
ions in the cascade region becomes comparable to the charge density on 
the anode wire, saturation occurs due to space charge effects and the 
strict proportionality ceases to exist. Because of the delta ray problem 
in high energy cosmic ray experiments it is desirable to operate MWPC's 
well within the proportional region and not at such a high gas gain that 
saturation occurs. In this respect it is important to note that the 
range of strict proportionality can be influenced to some extent by the 
choice of the MWPC gas^. 

The results of our MWPC investigation which are described in this report 
have all been obtained with "normal" gas mixtures as opposed to "magic" 
gas mixtures. MWPC's operated with normal gas show a region of strict 
proportionality before going into saturation and eventual breakdown. 
Saturation is due to space charge effects and breakdown is due to photon 
processes. With magic gas mixtures it is possible to prevent chamber 
breakdown and to obtain considerably higher gas gains than with normal 
gases, however, without maintaining proportionality. 

7 

R. Bouclier et al. who first reported measurements with magic gas ob- 
tained the gas by adding small amounts of highly electronegative freon - 
13B1 to their argon/isobutane mixtures. Extensive investigation of MWPC's, 

g 

operated with magic gas have been reported by J. Lacy. According to 
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Bouclier et al. it is the formation of heavy negative ions due to the 
electronegative nature of freon - 13B1 that seems to neutralize the pos- 
itive ions of the space charge, and thus makes it possible to achieve 
higher gas gains. Lacy believes and has indications that a basic change 
in the discharge process is responsible for the much higher gas gain. 
Indications are that only the first or first few electrons reaching the 
anode region are responsible for initiating the cascade process and that 
the high gas gain might be provided by a photon propogation mechanism 
which is controlled and confined by the photon absorbtion properties of 
the freon. Lacy has pointed out some attractive properties of MWPC's 
with coaxial delay-line readout operating at high gas gain with magic gas. 
At TRW we have concentrated on normal gases because extremely high gas 
gains are not necessary to achieve the desired spatial resolution, and 
because of some difficult problems that arise from magic gas operation. 

One severe problem arises from the loss of the ability to discriminate 
against delta-rays which is based on the proportionality of the chamber. 
Another potential problem is the fact that MWPC gases with organic con- 
stituents are subject to decomposition which can cause contamination of 
the MWPC and limit its lifetime. 

It is well known that the performance of MWPC's deteriorates due to con- 
9 10 

tamination 1 . Two forms of contamination exist: Contamination of 

the MWPC gas and contamination of MWPC wires. The MWPC gas can become 
contaminated with outgassing products from the materials that are used in 
manufacturing of the chambers. This requires the use of low outgassing 
materials and careful cleaning techniques. In the BCR-5 experiment the 
materials had to be limited to non-magnetic materials (except for the 
anode and cathode wires) . A way to overcome the gas contamination problem 
is to constantly flow fresh gas through the chambers (BCR-5 approach), 
or to dump the old gas and refill with new gas in periodic intervals. 

Another and potentially more troublesome problem is the possible contam- 
ination of MWPC wires. It has been observed that solid and liquid de- 
posits form on the wires after the accumulation of a large number of 
counts in gases containing organic constituents, especially hydrocarbons. 
The deposits apparently are polymerization and/or disintegration products 
of the organic constituents. The level of contamination is a function of 
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the total accumulated charge and thus depends on the number of counts and 
the gas gain. While the flux of high energy cosmic rays itself is not 
high enough even in two years to reach the critical number of counts, it 
has to be remembered that MWPC's are constantly sensitive so that the 
total number of counts is actually appreciable higher due to low energy 
background radiation. In this respect it is particularly important that 
the chambers can be turned off when the spacecraft passes through the 
South Atlantic Anomaly of the radiation belt with its high fluxes of low 
energy electrons and protons . 

The contamination problem not only requires careful material selection for 
the construction of MWPC's for space applications but also makes it de- 
sirable to operate the chambers with gases not containing organic con- 
stituents . 

3 . 1 Readout Concepts 

Several techniques have been used to read out the spatial information 
contained in a MWPC-signal after the passage of, for instance, a high 
energy cosmic ray particle or the conversion of an x-ray. The readout 
method which is actually chosen depends on several parameters, the most 
important of which is usually the spatial resolution that is required. 

Single wire readout systems and coding techniques are being used in large 
area MWPC systems if spatial resolution of a few millimeters is sufficient. 
This is often the case in large accelerator experiments and also in certain 
cosmic ray experiments like the ionization spectrometer experiment ACR-6 
on HEAO-A. In a single wire readout system each wire is connected to its 
own amplifier and subsequent chain of post amplifiers, gates, discrimin- 
ators and storage elements^. Readout electronics then determines the 
address of each wire that has responded. The uncertainty with which the 
position of a track can be determined in this way is equal to the anode 
wire spacing. In other words, a single wire readout system quantizes the 
position information and the quantum size corresponds to the anode wire 
spacing. An important advantage of this system is its capability to read 
out many simultaneous events. The disadvantage of such a system is the 
large number of electronic channels, with its impact on cost, power, weight, 
volume and reliability, especially if a large number of MWPC's has to be 
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used. For space applications a single wire readout system, therefore, 

becomes very easily prohibitive. For accelerator experiments the single 

12 

wire readout system is still feasible , even if it involves thousands 

of channels, because of the much more relaxed requirements on weight, 

power, volume and reliability which allow the cost to be kept at a few 

13 

dollars per channel 

The amount of electronics can be reduced if coding techniques can be 
used in which the signal coming off the anode wires are split and routed 
into several electronic channels which then can be decoded to give the 
address of the wire that has responded. The disadvantage of such coding 
systems is the reduced signal to noise ratio because of the large number 
of signal splits which could require operation of the MWPC at uncomfortably 
high gas gains. An additional disadvantage is the reduced capability to 
analyze several simultaneous events, and in large MWPC systems the amount 
of electronics will still be large. Single anode-wire readout systems 
become totally impractical if spatial resolution on the order of a 
fraction of a millimeter is required due to limitations in anode wire 
spacing. 

This problem can be overcome by making use of the induced charge on the 
cathode planes of a MWPC. The induced charge distribution peaks exactly 
above that point of the anode wire that collects the electron cascade, 
which makes it possible to retrieve position information perpendicular 
to the anode wire direction with considerable accuracy. Delay-line read- 
out systems have been developed which can interpolate the peak of the 
induced charge distribution to within .1 mm to .2 mm, i.e. considerably 
more accurately than the anode wire spacing. 

In the literature and also in this report we frequently use the term 
"spatial resolution" when we actually mean the accuracy of position deter- 
mination. Spatial resolution more accurately should describe the ability 
to resolve two closely spaced particle tracks or x-ray events. Spatial 
resolution defined in this way would, of course, depend on the accuracy 
of position determination, which is characterized by the width of the 
measured position distribution. 
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In a delay-line readout system the delay-line is coupled to either one or 
both of the cathode planes with the cathode wires running perpendicular 
to the anode wires. The position of the peak of the induced charge distri 
bution is determined by measuring the difference in arrival time of the 
delay-line signal at both ends of the line. Figure 3-3 shows the basic 
block diagram of a MWPC with delay line readout. 

It is obvious that such a system significantly reduces the amount of 
electronics necessary and that it can therefore reduce the complexity of 
a MWPC - system even when very good spatial resolution is not required. 

At TRW we have therefore, also designed and tested MWPC's with delay-line 
readout of the anode wires. The disadvantage of a delay-line readout 
system that replaces a single wire readout system is it's reduced capa- 
bility of handling a large number of simultaneous events. However, in 
many applications the handling of two or three simultaneous events is the 
most that is required, and this can be achieved with delay- line systems 
with some added complexity in the electronics. 

The spatial resolution is determined by the accuracy with which the peak 
of the induced charge distribution represents the position of the ionizing 
event in the MWPC gap, and by the timing accuracy of the delay-line read- 
out system. The electronic resolution or the jitter in timing At of the 
two delay-line signals is basically proportional to the risetime of the 
delay line pulse and to the signal to noise ratio of the output of the 
delay-line amplifiers: 

At 'v AR (RMS) 

As will be pointed out in the delay-line analysis in Section 6, the 
risetime is determined by the bandwidth of the line and/or the risetime 
of the MWPC current pulses. The signal level A s is a function of the 
chamber's gas gain and of the delay- line impedance. The impedance can 
also contribute significantly to the noise level with the other major 
contribution coming from the delay-line amplifier. 

14 

Grove et al. developed a distributed delay-line consisting of a solenoid 
al winding on a rectangular fiberglass core which forms the inductance. 

The distributed capacitance is formed by copper strips running along the 
core. To compensate for the inherent dispersion of the line the floating 
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FIGURE 3-3 Block diagram of a delay-line readout system. 

G = Delay-Line Amplifier, CFD = Constant Fraction 
Discriminator, THC = Time to Pulse Height Converter, 
PHA = Pulse-Height Analyzer, Z D = Delay-Line Impedance 



patch technique was used. Strips of Mylar with thin metallic strips on 
it was glued to the line. The delay-line was capacitively coupled to one 
of the cathode planes of the MWPC by clamping it to a flat PC-board with 
copper lands to which the cathode wires were connected. At TRW we have 
built and used a number of similar delay-lines. 

While excellent resolution (0.1 mm with collimated Fe-55 sources) and good 
linearity can be achieved with this delay-line, if it is carefully con- 
structed, it has a number of disadvantages which make it less desirable 
for large area MWPC systems in space applications. Considerable parastic 
capacitances in large MWPC's limit the coupling efficiency between chamber 
and delay-line to only 5% - 10%, which requires operation of the chamber 
of a much higher gain than would be necessary for a delay-line with the 
same impedance but larger coupling efficiency. There are also signifi- 
cant end effects of the line which are difficult to compensate and to 
control, and the adequate packaging of the line would result in a rela- 
tively heavy mechanical structure. Since the original work of Rindi et 
15 14 

al. and Grove et al. various other methods of implementing delay- 
line readout systems have been developed. 

A basically very simple way to implement a delay line with excellent 
coupling efficiency, linearity and no end effects is to connect adjacent 

cathode wires with pieces of coaxial cable. The use of coaxial delay- 

16 8 
lines has been reported by D. P. Boyd et al . of Stanford and by J. Lacy 

of Houston. The main advantage of coaxial delay-lines is their large 
bandwidth when used in connection with MWPC's operating with magic gas 
which results in considerably faster current pulses than those obtainable 
with normal gas for penetrating particles. Because of the large bandwidth 
of coaxial delay-lines the chamber risetime can be maintained and excel- 
lent spatial resolution can be achieved even with the reduced impedance 
of the line (50ft to 100ft as compared to close to 1 Kft for the line of 
Grove et al . ) 

There are again, however, severe disadvantages for space applications. Be- 
cause of the fast risetime (which allows smaller delays between wires than 
with the distributed line) fast timing electronics has to be used, which 
can lead to prohibitively large power consumption especially when a large 
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number of MWPC's has to be used. Another problem is the packaging of 
coaxial delay- lines, which become bulky and heavy in large systems. In 
order to improve the packaging, printed circuit delay-lines have been 
developed (D, P, Boyd et al^.) with characteristics similar to coaxial 
cable . 

A third method of implementing a delay-line readout system, which was 

17 18 

first mentioned by Charpak et al . , has been investigated by Lee et al. 
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and Gilland and Emming. In this case the chamber electrodes themselves 
form the delay-line (integral delay-line). In Lee's chamber the cathode 
consists of a bifilar helical winding of which one winding is grounded at 
every turn. The distance between windings was 0.5 mm and the cathode wire 
diameter was 50y. This resulted in such a high attenuation of the line 
that it was necessary to connect electrical tabs every 3.65 cm and to 
measure signals at these points. The obvious disadvantage for space ap- 
plications of the otherwise simple approach is the considerable complexity 
and large amount of readout electronics due to the many taps along the 
line for large area chambers. 

A somewhat different approach was taken by Gilland and Emming who intro- 
duced sense planes located symmetrically between the anode and cathode 
planes of a MWPC and which were wound as a continuous coil. The coil in- 
ductance and its distributed capacitance to the anode and cathode planes 
determine the delay line characteristics. Again, while the basic approach 
is elegant there are potential disadvantages of the readout system in space 
applications, especially for cosmic ray experiments. The major problem as 
we see it is the fact that the delay-line parameters and the achievable 
resolution are strong functions of the geometrical configuration of the 
chamber. This allows little flexibility in chamber design and can lead 
to difficult mechanical problems in chamber construction. This is 
especially true if the gap width of the chamber has to be kept small be- 
cause of potential problems due to delta-rays and inclined particles in 
cosmic ray experiments. 

At TRW we have developed a fourth approach to implement a delay-line read- 
out system that overcomes most of the difficulties of the other systems. 

Our delay-line which is called the "lumped element delay-line" consists of 
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individual L's and C's, which are directly coupled to both cathode planes 
and which determine the delay-line characteristics. A complete analysis 
of the delay-line is given in Section 6. 

In this analysis we first calculate the impulse response of a lumped ele- 
ment delay-line for a 50 cm chamber with the geometrical configuration of 
our standard chamber (Figure 3-1) operated with Ar/C02 gas, using the 
known design equations for the delay-line (e.g. equations for the risetime 
and delay per section, for the risetime and delay for n sections, and for 
the characteristics impedance) . Knowledge of the impulse response is 
necessary to calculate the delay-line output which is the weighted sum of 
the individual, induced cathode current pulses. The current pulses have 
been modeled according to the shape of observed current pulses for 5.9 
KeV x-rays from Fe-55 sources. When the delay line output pulse is known 
it is possible to calculate the spatial resolution (electronic resolution 
only) of the readout system in Figure 3-2, taking into account the rise- 
time and impedance of the delay-line, the amplitude of the delay-line 
signal, the noise of the preamplifiers, and the timing jitter of the CFD's. 
The analysis shows that an optimum interwire delay or delay per delay-line 
section exists for optimum spatial resolution of the particular MWPC under 
consideration. It is also shown considering the electronic system only, 
that a spatial resolution of at least lOOy for Fe-55 x-rays can be achieved 
at gas gains for which the chamber is in the proportional region. 

The advantages of the lumped element delay-line can be summarized as fol- 
lows : 

1. The delay-line characteristics are independent of the MWPC design 
(MWPC parasitic capacitances are taken into account in the delay- 
line optimization analysis). This allows flexibility in chamber de- 
sign and independent optimization of chamber and delay-line. 

2. Direct coupling of both cathodes into the delay line assures maximum 
coupling efficiency. 

3. Compact construction is possible, providing minimum weight and volume. 

4. Uniformity of the line makes it amenable to optimization studies. 
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3.2 Summary of Parameters Which Influence Spatial Resolution . 

It is possible to achieve good spatial resolution with large area MWPC's 
by using delay-line readout systems which measure the position of the 
peak of the induced charge distribution on the cathode planes of a MWPC. 

The accuracy with which the position can be determined depends on a number 
of parameters which determine how accurately the peak represents the 
position of the ionizing event and parameters which determine how ac- 
curately the electronic readout system can determine the position of the 
peak. These parameters are: 

1. Parameters associated with the gas multiplication process 

a. Diffusion of electrons and positive ions 

b. Spread of the discharge along the anode wire 

2. Parameters associated with the penetrating particle 

a. Inclination of the particle's track 

b. Delta rays 

3. The signal to noise ratio which, in a given system, depends on the 
gas gain. 

A detailed discussion of these parameters follows in Section 5 through 8 
of this report in conjunction with theoretical investigations and experi- 
mental results. The effect of electronic noise will be discussed in 
Section 6. 

In principle it is desirable to minimize electronic noise to the point 
where it is not the determining factor that limits the achievable spatial 
resolution, and our analysis will show that this is possible without having 
to operate the MWPC at such high gas gains (large signal to noise ratio) 
that one leaves the strictly proportional region. 

Attempts have been made to calculate the effect of electron diffusion and 
spread of the discharge along the anode wire on spatial resolution. Indi- 
cations are that both effects do not prevent a resolution of better than 
lOOy . Our measurements show, however, that some limiting effect does 
exist which results in a plateau in the measured resolution instead of a 
continuous improvement with improving signal to noise ratio. The con- 
clusion that has to be drawn from this apparent discrepancy is that it is 
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difficult to accurately calculate the effects of electron diffusion and 
discharge spreading because of some simplifying assumptions that have to 
be made and because some of the gas dependent parameters entering the cal- 
culations are not very well known. 

To estimate the effect of delta rays on spatial resolution we have carried 
out Monte Carlo calculations to simulate the effect of delta rays. The 
results of the computer calculations show that low energy delta rays 
broaden the main core of ionization which can definitely limit the achiev- 
able spatial resolution for high-Z nuclei (Z>4) to greater than the de- 
sired value of lOOy to 200y. For Z <_ 4 the delta ray effect is in the 
order of other limiting effects, and particularly for minimum ionizing, 
singly-charged particles (electrons and protons) delta rays are not import- 
ant . 

In cosmic ray experiments particles enter the spatial detectors on the 
average under some angle of inclination. The performance of MWPC's with 
collimated radioactive sources and normally incident particles . is , there- 
fore, not necessarily representative for the performance under the condi- 
tions of the actual experiment. We have measured the dependence of the 
spatial resolution on the angle of incidence and the result is that the 
resolution deteriorates with increasing angle of incidence. 

For measurements with collimated radioactive sources the collimation needs 
to be good enough so that the source width does not limit the achievable 
resolution. While good collimation of an Fe-55 source is no problem, it 
is more difficult to adequately collimate electron sources because of 
multiple scattering of the electrons. An arrangement to limit multiple 
scattering will be described in Section 7. 

4. GAS GAIN OF MWPC's 

In order to optimize MWPC configurations it is desirable to be able to cal 
culate the expected gas gain as a function of the MWPC geometrical configu 
ration, operating parameters and chamber gas. Several methods are avail- 
able from the literature to calculate the gas gain of single wire cylin- 

20 

drical proportional counters. We will show, however, that the method 
usually used fails to predict the gas gain of MWPC's of different con- 
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figurations (e.g. different anode wire diameters). 

We will then describe a straight forward calculation which is based on 
the knowledge of the first Townsend coefficient. The problem with this 
method is that the Townsend coefficient is not known for most gas mixtures 
used in MWPC's. We have, therefore, measured the Townsend coefficient 
for 80% Ar + 20% CO2, and we show that the calculations using the Townsend 
coefficient are self consistent and predict the gas gain for different 
geometrical configurations with fairly good accuracy. This points out 
the need for accurate measurements of the Townsend coefficient for the 
various MWPC gases over the E/p ranges encountered in MWPC's. 

4 . 1 Calculations Based on a Formula Given by Rose and Korff^ 

The calculation of the gain of a cylindrical proportional counter is 
usually based on a formula derived by Rose and Korff as given by Curran 
and Craggs: 

1/2 1/2 

In G = 2(2C V 0 a r x ) {(V Q /V S ) -1} (4-1) 

G = gas gain 

C = capacitance per unit length of anode wire 
a = constant 
r^ = anode wire radius 
V Q = applied voltage 

V s = threshold voltage, onset of gas amplification 

This formula can be applied without changes to a MWPC since the field 
around the anode wire is cylindrical and the gas multiplication process 
takes place in that region. 

Figure 4-1 shows the geometrical parameters of a MWPC: 



Cathode at GND 
Anode at V Q 


FIGURE 4-1 
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We define 


X = In 


f— ) 
sin h s 

u r'RT\ 

sin h ( — L ) 
v s 


or approximately 


X = It - In llILL 
s s 

The electric field around the anode wire is 


E (r) = 


Xr 


The charge per unit length of anode wire is 


* = 2x 


and the capacitance per unit length is 


C = 9 — 

V 0 2X 


(4-2) 


(4-3) 


(4-4) 


(4-5) 


With that we can write 

In G = 2(V 0 a ri /X) 1/2 {(V 0 /V s ) 1/2 -1} (4-6) 

which gives the gas gain as a function of the geometrical configuration 
and operating voltage of the MWPC. Equation (4-6) contains the constant 
"a" and the parameter V s which have to be determined. 

The constant "a" is gas dependent and is defined through 


0 


i 


1 

nX^ 


a(V-Vi) 


is the ionization cross section, X^ is the mean free path between 

3 

ionizing collision, n is the number of atoms per cm , is the ionization 
potential, and V is defined by the energy eV of the ionizing electron. 

The value of "a" thus expresses the rate of increase of the ionization 
cross section with energy. This has been measured for a number of gases 
but unfortunately not for the gas mixtures usually used in MWPC's. 
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At the threshold voltage V g the electric field at the anode wire is such 
that the ionization potential can be reached by the electron in a mean 
free path 



and 


V = Xr l V i 

s Ai 

V s is not only gas dependent but depends also on the geometry of the 
chamber. 

Because of insufficient knowledge and lack of measurements of gas para- 
meters for proportional counter gases both "a" and V g cannot be deter- 
mined accurately. V s also cannot be measured by pulse techniques because 
the signal level at the threshold voltage is too low to be recovered out 
of the noise. 

Good agreement between the predictions of equation (4-6) and experimental 
measurements has been achieved with suitable adjustment of the two para- 
meters "a" and V s for one set of conditions, i.e., for a particular gas 
and a given geometrical configuration. If the geometrical configuration 
is changed, e.g. by changing the anode wire diameter, equation (4-6) how- 
ever, does not predict the gas gain dependence with any reasonable accuracy 
To show this we used three MWPC's which were identical except for the anode 
wire diameter which were lOOy, 50y, and 20y and we measured the gas gain 
as a function of V Q for all three chambers in 80% Ar + 20% CO2. Figure 
4-2 shows the measured gas gain curves. If we fit equation (4-6) to the 
20y data and predict the lOOy data we get the results plotted in Figure 4-2 
As one can see, there is poor agreement between the calculated and measured 
curves. The discrepancy has to be due partly to the inaccurate knowledge 
of "a" and V s and partly to the simplifying assumptions in the derivation 
of equation (4-6) . It was assumed that the electron energy after each 
collision is always equal to eV^ and fluctuations in energy loss have been 
neglected. 
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GAS GAIN 
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FIGURE 4-2. Gas gain for 5.9 KeV x-rays in 80% Ar + 20% CO 2 
in chambers with different anode-wire diameter 
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We have attempted another approach to calculate the gas gain which 
parallels the one of Rose and Korff, but which uses a different set of 
simplifying assumptions. It leads again to an equation with two adjust- 
able parameters which allow a good fit to one set of measurements, but 
which also do not predict another set of measurements taken with a dif- 
ferent geometry. 

4 . 2 Calculation Based on Townsend Coefficient 

A straight forward way of calculating the gas gain of a proportional 
counter is possible by using the first Townsend coefficient, a, which 
characterizes the gas multiplication process. The Townsend coefficient 
is a function of the electric field E and depends on the gas mixture. 

When ot is known the gas gain can be directly calculated from the formula 

r2 

G = exp {/ a(E) dr} (4-7) 

ri 
or 

r 2 

In G = / a (E) dr (4-8) 

r l 

r^ is the anode wire radius and T 2 would be the cathode radius in case of 
a cylindrical counter. In a MWPC T 2 has to be large enough so that the 
multiplication process occurs inside a cylinder with radius T 2 around the 
anode wire. The Townsend coefficient has been measured as a function of 
E/p (p = gas pressure) for a number of pure gases and some gas mixtures, 
but again not for most of the gas mixtures used in MWPC's. 

The value of a can be derived from measured gain-voltage curves by dif- 
ferentiating equation (4-8) : 

r 2 

InG = f a(E) dr 
r l 

E - E(r) - £ 

InG = / Xr la(E) ^ ^2 (4-9) 

V X E 

v o 

xTJ 
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After differentiation: 


din} / dV 0 = If. ♦ 


«Cv^) • - "Cfe.) • P- 


v„ 


L Xr 


o 1 V o Ar 2 o 

As we will show below, because of the rapid decrease of a with E at low 


Xr-; 


values of E the term af^) vp- is less than 1% of the term a(r^) and 

, v o o 

can be neglected. 

The field at the anode wire is 


F = — o. 
Xr n 


and close to the cathodes it is 


E = 


V n TT 

Xs 


For our three chambers with r^ = 10 P> 25 P» and sou we had s = .4 cm, 
6.48 <_ X <_ 8.09, and 1500 V q <_ 3000. Near the anodes we get 

E/p > 105 V cm -1 Torr" 1 


Near the cathodes we get 


a/p > 1 

a/p • r^ > 50 x 10 


-4 


E/p < 4.8 V cm * Torr * 


-4 


a/p < 10 

a/p • L < 0.5 x 10 
which is less than 1% of the anode value. 

We can now express the Townsend coefficient as a function of the electric 
field at the anode wire as 


r\i /v -> dlnG InG 

a (V 0 /Xr 1 ) = _ - _ 

o o 


(4-10) 
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Before we calculate a/p as a function of E/p from our gas gain measure- 
ments with 20 p, 50y and lOOy anode wires, we have to remember that the 
measured values do not represent the total gas gain but the relative gas 
gain as pointed out in Section 3. Therefore, before we calculated a/p 
we calculated the absolute gas gain from the measured relative gas gain. 
The fraction of the total charge which is represented by the relative 
gain was given in Table III-l, which shows the time development of the 
MWPC pulse for the three different chambers. Figure 4-3 shows a/p as a 
function of E/p as derived from our adjusted gas gain measurements. 

To calculate the gas gain for 20y, 50y and lOOy diameter anode wire we 
used a linear approximation of a/p, which is the solid curve in Figure 
4-3. With that it is easy to integrate equation (4-8): 


Approximation for a: 


In G = / °a(r)dr 
r l 


a = D(E-E 0 ) for E > E Q 


a = 0 

Close to the anode we have 


for E < E 
— o 


E = V /Xr or, r = V /XE 
o o 

and r = V /XE which is the radius at which a = 0. 
o o o 

With that is follows: 

r o 

InG = f D(E(r) - E Q )dr 


V XE o 

= 0 f (V 0 /Xr - E 0 ) dr 


V„/XE. 

= D {V /X In r - E r) 
o o 

r l 
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FIGURE 4-3. Townsend Coefficient a/p v.s. E/p 
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which allows one to calculate the gas gain as a function of the geometri- 
cal configuration and operating voltage. 

The calculated gas gain curves using equation (4-11) for 20y, 50y and 
lOOy anode wires are also plotted in Figure 4-2 using E Q /p = 30 V cm ^ 

Torr E q (STP) = 22,800 V and D = .0262 V \ The plotted values are 
the relative gas gains which were obtained after adjusting the absolute 
gain of equation (4-11) according to the percentage numbers given in 
Table III-l. 

The agreement between measured and calculated curves is quite good for 
20y and 50y. The discrepancy for lOOy could indicate that the chosen 
approximation for a/p does not hold for low E/p values which is expected. 
However, the internal consistency of this method has been demonstrated 
and it would, therefore, be desirable to have accurate measurements of 
a/p available for the commonly used MWPC gases over the E/p range normally 
encountered in MWPC's. 

4.3 Gas Gain Measurements for Different Gases 

We have measured relative gas gain curves for four different gases in 
our standard MWPC which was operated in a flow-through mode. The 
measurements were made using a radioactive Fe-55 source. The anode 
wires in the center of the chamber were interconnected and then connected 
to an Ortec 109 PC charge sensitive preamplifier (see Figure 9-2 which 
shows the setup of our MWPC's for measurements in a heavy ion beam at 
the Bevatron) . The amplifier output signal was routed through an Ortec 
410 main amplifier into a pulse height analyser. The double differ- 
entiation, single integration time constants of the 410-amplifiers were 
set at 1 ysec (this applies to all our measurements in this report). 
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Four different gases have been used: 


1 . 

90% 

Ar 

+ 

10% 

co 2 

2. 

80% 

Ar 

+ 

20% 

co 2 

3. 

75% 

Ar 

+ 

25% 

Isobutane 

4. 

81% 

Ar 

+ 

14% 

Xe + 5% CO 


The measured gas gain curves are shown in Figure 4-4, in which we have 
plotted the number of electrons at the amplifier input v.s. the operating 
voltage of the chamber. The number of electrons represents the fraction 
of the total charge which is accumulated in 1 ysec. To calculate the 
relative gas gain one has to divide the number of electrons by the number 
of primary electrons in the chamber gap. This was done in Table IV-1 for 
the measurements in 80% Ar + 20% CO 2 . 

For comparison purposes we have also included some of the measurements of 
21 

Charpak et al . in Figure 4-5. Our measurements are consistent with the 
conclusions that can be drawn from Charpak' s more extensive results: 

1. The operating voltage necessary to achieve a certain gas gain in- 
creases with increasing percentage of quenching agent in the chamber 
gas. 

2. The onset of deviation from strict (i.e. linear) proportionality 
shifts to lower gas gains with increasing percentage of quenching 
agent . 

From an operational standpoint it might seem, therefore, desirable to 

use a chamber gas with a minimum amount of quenching agent. According 

21 

to results of Charpak et al . , which are shown in Figure 4-6, one has to 

be careful, however, in reducing the amount of quenching agent because at 
low percentages of quenching agent the voltages at the limit of the 
strictly proportional region and at the limit of chamber breakdown merge 
together. 

Our measurement with 81% Ar + 14% Xe + 5% CO 2 also confirms the results 
of R.Z. Fuzesy et al. that mixtures with Ar/Xe/C02 produce larger 
regions of strict proportionality than Ar/Q^ or Ar/ Isobutane mixtures. 

A result of practical importance is also the fact that Ar/Xe/C02 mixtures 
achieve high gas gains at considerably lower voltages than the other 
mixtures . 
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NUMBER OF ELECTRONS 



FIGURE 4-4. Gas gain for 5.9 KeV x-rays in different gases 
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TABLE IV-1 


GAS GAIN IN 80% Ar + 20% C0 2 


o Energy per ion pair: 27.2 eV 

o Fe-55 main peak: 5.9 KeV 

o No. of ion pairs in Icm 80% Ar + 20% CO 2 : 217 


HV No. of Electrons Gas Gain 


1400 

1.52 

X 

10 4 

7 

X 

10 1 

1500 

2.82 

X 

10 4 

1.3 

X 

10 2 

1600 

5.4 

X 

10 4 

2.49 

X 

10 2 

1700 

1.08 

X 

10 5 

4.98 

X 

10 2 

1800 

2.33 

X 

10 5 

1.07 

X 

10 3 

1900 

5.02 

X 

io 5 

2.31 

X 

10 3 

2000 

1.08 

X 

10 6 

4.98 

X 

10 3 

2100 

2.28 

X 

10 6 

1.05 

X 

10 4 

2200 

4.93 

X 

10 6 

2.27 

X 

10 4 

2300 

9.42 

X 

10 6 

4.34 

X 

10 4 

2400 

1.64 

X 

10 7 

7.56 

X 

10 4 

2500 

2.21 

X 

10 7 

1.02 

X 

10 5 
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FIGURE 4-6. Limit of proportionality v.s. limit of breakdown from 
Charpak's measurements (Ref. 21) 
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4.4 Space Charge Effects 


It has been pointed out before that the cloud of electrons and ions 
produced in the multiplication region near the anode constitutes a 
space charge that modifies the electric field configuration. This 
leads to a reduction in gain from the value that would be expected 
as a result of a calculation, such as we previously carried out, 
which neglects the presence of the space charge. The actual phenom- 
enon is a complex time and space dependent process that would be very 
difficult to analyze. We have performed some calculations on the 
space charge effect with very simplistic assumptions. The most im- 
portant of these are to completely neglect any time dependence, to 
use cylindrical geometry, and to use the simple linear expression for 
the first Townsend coefficient that was previously developed. As we 
shall see, the value of these calculations is to enable us to infer 
some limits on the effective width of the. charge distribution along 
the anode wire. Because the space charge effects depend on charge 
per unit length along the anode they are inversely proportional to the 
length over which a particular amount of charge is spread. 

Unless otherwise specified, we will be dealing with the following case: 

Chamber configuration: d = 2 x 10 _ 3 cm 

L = 0.5 cm 
s = 0.4 cm 

Gas: 80% Argon - 20% CO2 

Source: 5.9 KeV Fe^ x-ray 

(i.e. - 220 electron-ion pairs primary 
ionization) 

We shall first calculate the reduction in gain near the onset of space 
charge effects. For our operating conditions this occurs at about 2300 
volts (See Figure 4-4). In the case where the effects of space charge 
are small, we calculate the charge distribution without taking account 
of space charge, and then estimate the effect on the gain of this distri- 
bution. If the calculated effect is small, the method is consistent; 
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if not, we shall have overestimated the reduction due to space charge. 

In the next section we will estimate the upper limit on the charge that 
can be delivered under the extreme condition of space charge saturation. 

At a radius, r, the number of electrons inside r is equal to the initial 
number electrons times the gain due to multiplication up to that radius, 
i.e. , 

r 

n_(r) = n Q exp{/ adr'} 
r o 

There is obviously also an equal number of ions outside of r. If the 
charge distribution has an effective length along the anode, w, the 
total charge per unit length inside r is, 

q' = q - -2— exp{/ adr'} 
w r o 

where, 

q = V q /2X 

is the original charge per unit length on the anode wire. The reduction 
in the electric field at r is consequently, 

2n e r 

AE(r) = 2Aq(r)/r = - — exp{/ adr'} 

W r D 

We will use the following simple form for the first Townsend coefficient, 

a = D (E - E q ) , with, 

D = 0.0262 V” 1 
E q = 22,800 V/cm 

For V Q = 2300 V, the radius at which a goes to zero is, 

r Q = V 0 /E q X = .0125 cm, and 
x 

AE(r) = - exp {/ D(V 0 /X) ( p - p-) dr'} 

r© x o 


2n^e 
w 


exp (In ^ - 1) } 

r X r r„ 


Using the numberical values above, we get. 


n n 0046 596r 

AE(r) = Ja (2.88 x 10~ 7 ){ (■ ■ - ■ - ) e } V/cm 
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The change in gain is given by 


AlnG = / °Aa(r)dr 
r l 


o 

= f DAE (r) dr 
r l 

Using the expression for AE(r), above, and plugging in the numbers, we 
get; 

AlnG = 1.8 x 10 ^ n Q /w 

In order to estimate w, we will use the measured relative gains for the 
main Fe^ peak (n Q = 220) and the escape peak (n Q = 110) . At 2300 V we 
observed a gain ratio of 1.2 therefore. 


ln(l .2) 


AlnG 

mam 


AlnG 


escape 


0.175 


1.8 x 10 
w 


2 


.02 

w 


This result can be used to establish an upper limit on w, because the 
actual time dependence of the process can only decrease the space charge 
effects. Hence, we have, 

w <_ 0 . 1 1 cm 

We can compare this effective width with what we would expect from • 
electron diffusion that was estimated earlier. The effective width of 
the lateral distribution is. 


w = m a = 0.08 cm 

The comparison with measurement above essentially amounts to examining 
the change in the space charge effects due to a change in the initial 
charge with the effective width assumed to be constant. Alternatively, 
we can make use of the observed variation of gain with angle of incidence 
of a penetrating particle for a chamber that is being operated in condi- 
tions for which space charge effects are significant. We have such 
measurements with C + ^ ions at the Bevatron with a chamber using 75% Argon 
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- 25% Isobutane at 0, 10, and 20 degrees incidence. All we need from 
the analysis above is the result that: 

AlnG « n Q /w 

In this case, in addition to the spread of the discharge along the anode, 
we have a known geometrical contribution to the effective width of the 
charge distribution which is equal to 2LtanO. The number of primary 
electron-ion pairs is proportional to secO due to the increased path 
length of an inclined track. The results are shown in Figure 4-7, in 
which we have plotted AlnG*cos 0 v.s. w The relationship should be 
linear. In addition to the measured points at 10° and 20°, we know that 
AlnG must go to zero for large 0 because the charge would be spread over 
a very large length. If we extrapolate to the value of AlnG observed 
at normal incidence we get 


.050 <_ w <_ .070 cm 

Again, we have no indication of an effective width of the discharge that 
is not consistent with that expected from electron diffusion. 

4.5 Space Charge Saturation 

As positive charge accumulates around the wire, the field near the wire 
decreases, and the field at larger radii increases, in such a way as to 
keep the integral from anode to cathode equal to the applied voltage. 

The limiting value to which the field strength near the wire can be re- 
duced is that value at which a decreases to zero and multiplication 
ceases. For the 80% At - 20% CC >2 mixture with our simple form for a 
this happens at E Q /p=30 V cm'^Torr 

The maximum space charge per cm that can be accumulated corresponds to a 
field which is equal to E Q out to some radius r Q , and then drops off: 

There can be no further multiplication beyond this point. There is an ac- 
cumulation of charge out to r , but none beyond because the field there 
has always been less than E Q . 

In order to simplify the calculation, we now approximate the MWPC by a 
cylindrical counter with the same anode wire radius, r^, and the same 
capacitance per cm, C. For our MWPC configuration we have. 
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C = 1/2 X. 


For a cylindrical counter, 

C = 1/2 In 

where r 2 is the cathode radius. Therefore, 

lnr 2 = X + lnr r (4-12) 


The value of r Q is determined by the condition that the integral of the 
field from r^ to r 2 be equal to the applied voltage. At the wire sur- 
face : 


= ll 


r 


1 


In the same way, the field at the surface of any cylinder of radius r is 
determined by the charge inside the cylinder: 

Tn 


f ^2irr' Q(r' )dr' = j E Q r 


where Q is the charge density. 

g 

The solution is, Q(r) = 


The space charge per cm is , 


and the total charge is 
The applied voltage is 

V 

o 

V 

o 


q = / 2m:' Q(r')dr' 

r-, 


= 2 E o ( V r l> 


q+q' = T E r 
n n 2 o o 


E (r -r ) + / 2 E r 
o o 1 t o o r 


E (r -r,) + E r In ^2- 
o o 1 o o r Q 


(4-13) 
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Substituting lnr2 from equation (4-12), we get, 


V = E (r -r,) + E r (X + lnr. - lnr ), or 
o o^ o o o v 1 o J ' 


r /r 
o 


1 


1 ♦ , 

1 + X- lnr /r, 
o 1 


(4-14) 


This equation can be solved for a given set of operating conditions (i.e. 
X, V q , r , and E q ) to obtain r Q . The maximum charge per unit length, q', 
can then be determined from equation (4-14). We have calculated q' as 
a function of V 0 for our chamber configuration with 20 micron anode wire 
diameter. The results are shown in Figure 4-8. The value of q' is 
seen to be comparable to the charge per unit length on the anode as 
would be expected. We can arrive at an estimate of the maximum gain 
realizable as a function of voltage under the assumption that the ef- 
fective width of the discharge is equal to that due to electron dif- 
fusion (= 0.08 cm) and constant, since: 


q' = n e G /w, or for Fe 
n o max * 


55 


X-rays, 


G =wq'/ne= .8x 10^ q' (esu/cm) . 
max n o n 

The result is shown in Figure 4-2 along with the measured gas gain. If 

the measured gas gain exceeded G we would have a clear indication that 
6 6 max 

the discharge must be spreading along the anode. Since, in fact, this is 
not the case we can only obtain an estimated lower limit on the effective 
width of the discharge by assuming that the maximum output signal that 
we measured corresponds to the space charge saturated gain. At 2500 V, 
we obtained an estimated total charge output of about 5 x 10^ electrons. 
Therefore, 


> (5 x IQ 7 ) (4 , 8 x 10~ 10 ) 
.6 esu/cm 


esu 


> 0.04 cm 


In conclusion, we have seen that the highly simplified analysis of space 
charge effects has placed the following limits on the effective width of 
the charge distribution of: 
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FIGURE 4-8. 


Charge per unit length and maximum charge per unit 
length v.s. high voltage V Q 
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0.04 cm <_w <_0.11 cm, 

for operating conditions above the linear region but probably well below 
complete space charge saturation. This width is roughly the same as the 
effective width due to electron diffusion in the drift region. There- 
fore, we see no indication of significant spreading of the discharge a- 
long the anode under these conditions. It is possible that in the com- 
plex real situation of the space charge phenomenon, the shape of the 
charge distribution or the time dependence of the MWPC current pulse is 
being distorted in a manner that causes a degradation of spatial resolu- 
tion. 

5 . SPATIAL RESOLUTION 

Position determination using MWPC's with delay-line readout systems is 
usually based upon measuring the position of the peak of the induced 
distribution on the cathode planes. The accuracy of the position 
measurement is conceivably influenced by a number of parameters, which 
were mentioned in Section 3.2. These parameters are: 

a. Diffusion of electrons and positive ions 

b. Spread of the discharge along the anode wire 

c. Inclination of the particle track 

d. Delta rays 

e. Signal to noise ratio 

5.1 Diffusion of Electrons and Positive Ions 

It turns out that the lateral diffusion of electrons and positive ions 
in the cascade region in the immediate vicinity of the anode wire does 
not contribute to errors in the position measurement because of the 
large number of these charge carriers in this region and because of the 
small radial extension of this region. However, the diffusion of the 
primary electrons in the drift region from the cathode to the anode 
until the cascade process starts can contribute significantly towards 
a position uncertainty. 

The r.m.s. diffusion distance (one-dimensional) is given by. 


<x> 
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where D is the diffusion constant and t is the drift time through the 
anode-cathode gap. The problem would be straightforward if measured 
diffusion coefficients were available. Again we find that unfortunately 
such data does not exist for the commonly used MWPC gases. The diffusion 
coefficient is known for pure Argon, but it is not possible to enter a 
value for MWPC mixtures because a principle effect of the quenching ad- 
ditive (e.g. CC^) is to radically reduce the mean electron kinetic 
energy which directly influences the diffusion coefficient. It is pos- 
sible to calculate the diffusion coefficient if the mean electron kinetic 
• 22 

energy is known. Hough indirectly arrives at an estimate of mean 
electron kinetic energy and diffusion coefficients for Argon-methane 
mixtures. For 90% Argon - 10% methane he estimates and experimentally 
verifies a FWHM single-electron resolution of approximately 1 mm. This 
value includes electronic noise and source width contributions of 0.3 mm 
each. The diffusion contribution that we infer from these results is 
0.65 mm FWHM or about 0.3 mm rms. We have estimated the mean electron 
kinetic energy in 80% Argon - 20% C0 2 to be about 3eV for the drift 
region electric fields in our normal operating conditions. This leads 
to an estimated rms diffusion distance of 0.24 mm for a 5 mm drift dis- 
tance . 

For n primary electrons the r.m.s. diffusion is 

a = <x>//n~ 

For a minimum ionizing, singly changed particle which produces about 100 
primary electrons in our chamber, the position uncertainty of the peak 
of the induced charge distribution due to diffusion thus is approximately 


a = 25y 


at normal incidence. 


The position uncertainty of 5.9 KeV Fe-55 x-rays, which produce about 220 
primary electrons, is according to our calculations a = i6y. This does 
not take into account, however, the lateral extent of the primary ioni- 
zation due to the conversion process of the x-rays into electrons. 
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The x-rays are almost entirely absorbed by the argon atoms. The proba- 
bility that absorption leads to the emission of an electron from the K- 
shell is: 


and for the L-shell: 


P K = 0.90 


The fluorescence yields are: 


P. = 0.10 

JL 


“K = 0-H 

= 0.01 


The probability of escape of the K x-ray is: 

£ K = °* 9 

Eighty-one percent of the time we get K-shell absorption followed by de- 
excitation via the Auger effect. This will result in the production of 
two 3-KeV electrons at the point of interaction. 

Eight-to-nine percent of the time there is one 3-KeV electron (escape 
peak) . 

Zero-to-one percent of this time there are two ~ 3-KeV electrons at dif- 
ferent points in the chamber (reabsorption of the Ar K x-ray) . The mean 
free path of the Ar K x-ray is X = 3.8 cm; the two points will, therefore, 
in general be several millimeters apart. 

Ten percent of the time there are one 5.65-KeV electron and one 0.24-KeV 
electron (L-shell absorption) . 

The range of a 3-KeV electron in the gas is about 0.1 mm. For the two- 
electron case, if the directions of emission are not correlated, the 
r.m.s. distance along one axis between the point of origin and the cen- 
troid of the ionization distribution is 20 microns. If there is only 
one 3-KeV electron, this distance becomes 30 microns. These two cases 
between them cover 90% of the events. Since the value of a, the r.m.s. 
position error, is in general determined on the basis of the width of 
the curve at the 50% or 60% point, it will not be affected by the other 
10% of the cases, it appears that electron ranges contribute about 20-25 
microns to the measured value of o for Fe^ 9 x-rays. 
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5.2 Spread of the Discharge Along the Anode Wire 


The symmetry of the induced charges on both cathode planes and in the 
neighboring anode wires indicates that the discharge spreads around the 
whole anode wire regardless of the position of the primary ionization. 
Therefore, with an anode wire diameter of 20 microns, the discharge spreads 
20 to 40 microns in each direction around the wire. It seems reasonable 
to assume that the discharge spreads at least this far along the wire. 

As a matter of fact, due to the diffusion of the primary electrons, 
which is <x> = 240p for a single electron, the spread is several times 
this distance. As long, however, as the discharge spreads equally in both 
directions, the peak of the induced charge distribution accurately re- 
presents the position of the primary ionization. This is true as long 
as the peak remains well defined and doesn't spread significantly. 

We can get a rough idea of the upper limit on the spread of the discharge 
along the anode by observing the spatial distribution of the induced 
signals on the cathode wires. The expected induced charge distribution 
varies with position as 


1 


where L is the distance from anode to cathode. The f.w.h.m. of this 
distribution is 2L or about 10 mm in our chamber. We have measured the 
induced charge distribution in the strictly proportional region of our 
chamber and find a f.w.h.m. of 11 mm. This difference can be accounted 
for by the contribution of electron diffusion discussed above. This 
indicates that the discharge does not significantly spread along the 
anode wire in the proportional region, and that it might not contribute 
significantly to uncertainties in position measurements. There is the 
possibility that at sufficiently high gains, space charge effects might 
significantly affect the spread of the discharge along the anode. We 
examined this possibility in more detail in Section 4.5 and came to the 
conclusion that the width of the charge distribution at high gas gain is 
still roughly the same as the effective width due to electron diffusion. 

The results of our spatial resolution measurements, which will be pre- 
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sented below, show, however, that the measured spatial resolution satur- 
ates with increasing gas gain or signal to noise ratio, which can be 
interpreted as an indication that our simplified analysis does not re- 
present the complex real situation at high gas gains, and that the shape 
of the charge distribution is distorted in such a manner as to degrade 
the spatial resolution. 

5.3 Inclined Trajectories 

If a particle enters a MWPC with an inclined trajectory (inclined in the 
anode wire plane) the discharge will spread out along the anode wire and 
the amount of spread will depend on the degree of inclination. Measure- 
ments we have carried out with collimated electron sources show that a 
serious deterioration of the spatial resolution occurs for inclined 
tracks. This does present a problem for cosmic ray experiments which 
accept inclined particles over a wide range of angles. One way to re- 
duce the effect of inclined tracks in MWPC's operated with normal gas 
is to reduce the gap width in order to reduce the length of the dis- 
charge along the anode wire. 

5.4 Delta Rays 

A limitation on the spatial resolution of MWPC's is the spatial distri- 
bution of the primary ionization produced by the passage of a charged 
particle through the MWPC gas. Specifically, secondary electrons (delta 
rays) are produced which have energies such that their range is greater 
than the magnitude of the uncertainties in position location that we have 
discussed above. On the average, delta rays are produced symmetrically 
around the particle's trajectory and should have no first order effect 
on the position of the peak of the induced cathode charge. However, the 
frequency of high energy delta ray production is low enough that the 
statistical fluctuations in delta ray production can lead to asymmetries 
in the induced charge distribution and thus to significant uncertainties 
in the position measurement of individual events. In addition, the prob- 
lem can be aggrevated if the MWPC is operated above the strictly propor- 
tional region, in which case the primary ionization does not multiply 
linearily. In effect, the ionization produced by delta rays outside the 
main core of ionization will be weighted more heavily. 
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In order to make a quantitative estimate of the potential severity of 
the delta ray problem, a computer program was generated that uses Monte 
Carlo techniques to simulate delta ray effects (see Section 8) . The re- 
sults of these calculations, which have been carried out for particles 
with an energy of 2.1 GeV/nucleon and different Z values, show that delta 
rays can have a serious effect on position measurements. Indications 
are that the irreducible delta ray effects are essentially of the order 
of or less than the spatial resolution we want to achieve (100-200 microns) 
for relativistic nuclei with Z < 4, and increases to about 1 mm for Z = 

26. 

5.5 Signal to Noise Ratio 

The electronic noise of the readout system and the resolution limitations 
due to the signal to noise ratio will be discussed in the next section. 

The analysis of the delay-line readout system shows that the signal to 
noise ratio at moderately high gas gains in the proportional region is 
not a limiting factor. 

5.6 Summary 

In summarizing the above discussion we can conclude that none of the dis- 
cussed effects seems to significantly limit the achievable spatial re- 
solution for minimum ionizing, singly charged particles at normal inci- 
dence and for Fe-55 x-rays. A resolution of u < 100y should be achievable 
in the proportional region of our chamber. We have to remember, however, 
the assumption that the spread of the discharge along the anode wire is 
insignificant in that region. 

Accordingly we can calculate the theoretical limit of spatial resolution 
for our chamber for 5.9 KeV x-rays for the conditions of our measurements: 

Diffusion: 16y 

Lateral extension of primary ionization: 20y 

Source width: 45y (see Section 7.1) 

Electronic Resolution: 45y (see Section 6.4 and 

Figure 6-10 for 2100 volts) 
Theoretical Resolution: o = /16^ +20^+45^+45^ = 70y 

The real limitations in spatial resolution in high energy cosmic ray 
experiments are posed by the delta-ray problem and by inclined trajectories. 
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6.0 LUMPED ELEMENT DELAY-LINE ANALYSIS 


6 . 1 Introduction 

Various MWPC readout systems and their advantages and disadvantages have 
been discussed in Section 3.1. It was pointed out that accurate position 
determination of an ionizing event can be achieved with delay line read- 
out systems which determine the position of the induced charge distri- 
bution on the cathode planes of the MWPC. 

Since the achievable spatial resolution is proportional to the signal to 

noise ratio of the readout system, delay lines which are directly coupled 

to both cathode planes are preferred over delay lines which, for example, 

are capacitively coupled to only one cathode plane. The integral delay 

19 

line of Gilland and Emming could be considered the ultimate realization 
of a directly coupled delay line concept since it uses the MWPC itself 
as the line. However, as we have pointed out before, the strong depend- 
ence of the characteristics and performance of the integral delay line 
upon the MWPC geometrical configuration leads to a significant loss in 
chamber design flexibility. In fact, the dependence upon chamber configu- 
ration could lead to constraints in chamber design which are unacceptable 
to the experiment's objectives. Another disadvantage of the integral de- 
lay line is the added mechanical complexity due to the introduction of 
sense planes between the anode and cathodes, and the requirement that 
relatively large diameter wires have to be used for the sense planes in 
order to reduce the line attenuation to reasonable numbers. 

The lumped delay line readout system presented here allows near perfect 
charge coupling efficiency, independent optimization of chamber and delay 
line parameters and spatial resolutions of the order of 100 pm (for Fe 
X-rays) with the chamber operating at a reasonable gain in the proportional 
region. The use of lumped element delay lines for MWPC readout is not new, 
having been proposed by Charpak, et al.^; however, an analysis has not 
been presented that identifies the key design parameters and predicts the 
achievable performance. This section describes an approximate analysis 
that may be used to choose optimum delay line parameters and to evaluate 
the effect of electronics noise on the achievable spatial resolution. 
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6.2 Delay-Line Analysis 


The creation of an avalanche in the small region surrounding an anode 
wire causes charge to be induced on adjacent anode wires and on the 
cathode planes. If we defined as the total charge on the active 
anode wire and Qj, the total charge induced on the neighboring anode 
wires, then the total charge induced on the cathode wires, Q, is given 
by 

Q = Q a - Qj (6-1) 

Qj is dependent upon the geometrical configuration of the chamber and 
can typically be 10 to 20% of Q^. The induced cathode charge may be 
expressed as 


Qn 

Q 


1 

TT 


L/S 


& 


+ N 


( 6 - 2 ) 


for wire cathode planes. In Equation (6-2), L is the anode plane to 
cathode plane spacing, S is the cathode wire spacing and N is the wire 
number (with N = 0 being that cathode wire directly below the anode 
charge) . 


The analysis of the lumped element delay line readout system requires 
a mathematical representation of the cathode current pulses. Measure- 
ments performed on a small chamber indicate that the anode current pulse 
(and hence the induced cathode current pulse) may be modeled as 


i (t) = Q ±2 e' t/T (6-3) 

For Fe^ X-rays in a 90% A, 10% CC^ gas mixture we measured a model pulse 
time constant, t, of 20 ns. The model pulse used in this analysis ex- 
hibits a risetime of approximately 11 ns, a FWHM of 49 ns and a peaking 
time of 20ns. The model pulse time constant for X-rays is actually not 
very sensitive to changes in.operating conditions (gas and high voltage) 
over our operating range, due to the concentration of the primary ioni- 
zation in a small, localized region. Non X-ray events can be expected 
to have different x's. 

The amplitude of the model cathode current pulse of Equation 6-3 is de- 
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pendent upon the location of the cathode wire with respect to the anode 

til 

charge (Equation 6-2). Therefore, the current pulse on the N cathode 
wire can be written as 

V*’ " 72 e ' t/T 


Q L/S t -t/r 

* (L/S) 2 + N 2 t 2 6 


(6-4) 


Consider a tapped delay line of infinite bandwidth and a constant group 
delay of t^ seconds between taps. Each tap is connected to its cor- 
responding cathode wire. The delay line output may be written as 


M 

i t‘) = £ % (t-Nt dl )e- (t - Nt dl)\ (6-5) 

N=-M T dl J 

where M wires either side of the zeroeth are assumed to contribute 
significantly to the output. The function U(t) is the unit step function 
defined as 

r 0 t < 0 

U(t) =i 

L i t >_ o+ 

(The delay from N = 0 to the end of the line has been ignored here since 
all wires are affected equally.) 


The effect of a bandwidth-limited, lossless delay line is to reduce the 
amplitude of the pulse and to spread it in time. As will be demonstrated, 
the bandwidth of the line as measured from an input applied to one of the 
taps is a function of the number of delay line sections between the tap 
and the line output. One may model the physically realizable delay line 
of Figure 6-1 (a) as shown in Figure 6-1 (b) . Here the individual in- 
puts are applied to ideal delay lines whose delay is a function of the 
distance between the input and the line output. Each delay line output 
is then applied to a filter representing the bandwidth of the line as 
seen by an input applied to the particular tap in question. The delay 
line output is therefore the weighted sum of the individual cathode 
current pulses. Knowledge cf the impulse response of the component 
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filters shown in Figure 6-1 (b) is required in order to calculate the 
delay line output. The following paragraphs will present a discussion 
of a lumped element delay line and an approximation to its impulse re- 
sponse . 

A lumped element approximation to a delay line having equal delay between 
taps may be realized by cascading constant-k T sections, as shown in 

Figure 6-2. The end termination networds are m-derived half-sections 
including to improve the matching between the line and R Q . 

23 

The delay line design equations are: 

t j = 1.13/LCf, risetime per section 

t^ = 1.07/LCr, delay per section 

t, = n t.., delay for n sections 
a al 

t = 3 /rT t , , risetime for n sections 
r rl 

R o = characteristic impedance 

Therefore, for fixed n the delay to risetime ratio is constant. (Note: 
Should a faster risetime be desired for the same delay, additional T 
sections may be added to the line between tap locations.) 

The theoretical impulse response of an equal element lumped delay line 
involves Bessel functions of high order. Convolution of the model pulse 
and the line's impulse response yields the line's output signal. This 
problem in itself is a time consuming computer operation and when one 
considers that it must be done for each cathode wire containing a "signifi 
cant" fraction of the charge (because of the dependence of the order of 
the Bessel function on wire number), it becomes apparent that obtaining 
the theoretical solution is very expensive. 

In light of the computational difficulties just described, the following 
simplified analysis will be employed. The impulse response of the delay 
line can be approximated in the form shown below (neglecting the fixed 
delay) : 
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The 10 to 90% risetime of the waveform of Equation (6-7) is calculated 
to be t r = 0.4822T. As before, the risetime of the line is given by 


t 

r 


Ft 


rl 


or 


3 

T = 2.19 /n t,. 

dl 


Convolving the model pulse, 


( 6 - 8 ) 


i(t) 


Qot -t/i 
2 e 


with the impulse response approximation, Equation (6-6) yields, after a 
considerable amount of arithmetic, the delay line output for a single in- 
put : 


[l- (l+~) e 


-t/r 


,e-i 


e' t/T ,t 8 2 -10 


i'(t)=< 0, 


t<0 


* ( ' U6-!) 2 ^ (i - j^f)] ,0<t<T 


(6-9) 


12a e”^ T -l-— ♦ 2 [1 - 5lA - ( 

T V t x 1+6^ Lt r2+i 


6 2 -i 


rt-T 


B 2 + l 


e 2 -i. t/ti, 

g7TT) e 


where u = — , 8 = 


and T is given by Equation 6-8. 


Let the line have Nj sections in total with the anode point charge lo- 
cated above the ^ th tap, as shown in Figure 6-3. 
o — 
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FIGURE 6-3 
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It will be assumed that the Ml + M2 wires located around the center of 
mirrored charge distribution contribute significantly to the output 
pulse. (Note that as long as Ml ^ - N 0 and M2 ^ N Q , we will have 

Ml = M2 = M) . The delay line output, i-p(t), may now be written (from 
Figure 6-1) as 


Ml 

i T (t) = E i' N (t-Nt )U(t-Nt dl ) 
N=-M2 


( 6 - 10 ) 


Substituting Equation (6-9) into 10 yields 


Ml 

,(t)= E 
-M2 


( ij. (t -Nt dl ) [U (t -Nt dl ) -U (t -Nt dl -T)}+i 2 (t-Nt dl )U (t -Nt dl -T)} 

( 6 - 11 ) 


where is given by Equation 6-2 and ij and i 2 are defined below: 


i 1 (t)=l ~(1+ ^-)e' t//x + coswt-Bsinwt) 

-tlx 


+ l+e^” 


c 1 e2_1 ) 

( x " B 2 +l J 


i 2 (t) 


■t/x. 


. r,, t-T, T/t , t 1 rt B~-l ,t-T B -1^ T/t-i i 

e {(1 + — )e -1- 7 + nfrt- p2TT-C— - FTT )e 


It should be noted that the output waveform described by Equation 6-11 
has been shifted in time so that the cathode pulse applied at section 
N = 0 occurs at t = 0 at the output. The effect of the bandwidth vari- 
ation with input position is treated by noting that 

T = 2.19 /N 0 + N t dl (6-12) 

The delay line, when electrically connected to the chamber, exhibits an 
increase in risetime over the isolated line. Computer simulation of the 
delay line indicates that intertap parasitic capacitance added by the 
printed circuit pattern interconnecting the chamber and line is respon- 
sible for this effect. In order to account for the reduced delay line 
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bandwidth, the value for T given by Equation (6-12) is increased. For 
our typical laboratory delay line (L = 4.7 yh, C = 20 pf, Ci nter t a p = 

2 pf) , the observed risetime is approximately twice that predicted by 

3 

the design equations and hence we will see T = 4.38 /Nq + N t^ when 
referring to the "loaded" case. Accompanying the increased risetime is 
an increase in the overshoot in the step response; from approximately 6% 
for the isolated line to 30% for the loaded line. Because of the 
smoothing effect of the multiple cathode inputs, the output waveform 
appears relatively clean. In some applications, however, the ringing 
may be objectionable and a suitable phase correction filter should be 
incorporated at the delay line output. 

6.3 Implementation of the Lumped Element Delay-Line 

The method of interconnecting the cathode planes and the technique of 
coupling them to the delay line warrants some comment. For illustra- 
tion, the chamber will be assumed to be 50 x 50 cm with a cathode wire 
spacing of 2 mm. The resultant 500 cathode wires must be connected to 
the delay line. It has been previously mentioned that in order to achieve 
the highest charge coupling efficiency, both cathode planes must be dir- 
ectly connected to the delay line. 

Figure 6-4 indicates the evolution of the adopted configuration. Figure 
6-4 (a) shows a straight-forward mechanization of the constant-k delay 
line which is further simplified to the configuration of Figure 6-4 (b) . 
Employing the self-inductance of the cathode wires and the anode-cathode 
capacitance (augmented with external capacitance, if necessary) yields 
the delay line-chamber system of Figure 6-4 (c) . Notice that in all 
three drawings the winding sense of the cathode wires, i.e. essentially 
bifilar, is such as to minimize the mutal inductance thereby maintaining 
a constant-k delay line. An m-derived delay line may be easily developed 
from the last figure as shown in Figure 6-4 (d) . Although the m-derived 
form (similar to that achieved in the integral delay line chamber) offers 
a small potential improvement in group delay constancy with frequency 
over the constant-k line, it does require that the chamber be so sized as 
to yield the proper amounts of self and mutual inductances. The constraint 
imposed here may require a compromise between chamber parameters (such 
as gap width) and delay line parameters. 
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Each of the chamber-delay line systems of Figure 6-4 have the cathode 
wires in series with the delay line elements. If one uses the smallest 
cathode wire size compatible with mechanical and electrostatic stability, 
then a serious attenuation and noise source results from the wire re- 
sistance. For example, 50 y diameter Cu-Be cathode wires (spaced 2 mm 
in the 50 x 50 cm chamber) contribute over 6000 ohms in series with the 
line. Of course larger wire sizes may be used at the cost of increased 
mass density across the MWPC sensitive area and possibly larger frame- 
width if additional tensioning of the wires becomes necessary. 

The adopted configuration of Figure 6-5 avoids the aforementioned problems 
while still allowing for direct coupling to both cathode planes. The 
cathode wire-delay line connection alternates between the upper and lower 
cathode planes in order to achieve a bifilar winding. Connection between 
the cathode wires and the discrete delay line components are made by a 
printed circuit board. 

6.4 Effect of Electronic Noise on Spatial Resolution 

The effect of electronics noise upon system spatial resolution is readily 
determined. For the purposes of this report the readout system is as 
shown in Figure 6-6 and consists of electronics viewing each end of the 
delay line. Each delay line output is amplified by a low noise voltage 
preamplifier. The preamp outputs are then fed to two identical constant 
fraction discriminators (CFD's). The individual CFD outputs then become 
the start and stop inputs to a pulse height to time converter whose out- 
put is then accumulated in a multichannel analyzer. 

The preamplifier is modeled as an ideal (noiseless) amplifier with current 

and voltage noise represented explicitly by the two noise sources (assumed 

2 2 

to be uncorrelated) of noise spectral densities i^ and e^ , respectively. 
Figure 6-7. The delay line and its termination resistance, Rq, has as- 
sociated with it t he J ohnson noise voltage represented by the noise 

2 

spectral density, i . The signal, derived from the superposition of the 
individual cathode pulses, is labeled i^,. 

The total equivalent simplifier output noise spectral power density is 
found to be 
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FIGURE 6-6. Experimental setup for resolution measurements 
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(6-13) 


2 = r .2 


■'no 


= G z {e A 


. 2 R n .2, ,,,2.,, , 
1A -g” + X R } ( V / Hz ^ 


where 

i R 2 = 2KTR Q and 


where K is Boltzman's constant and T is the absolute temperature (°K). 

The bandwidth required to pass the delay line output signal is a function 
of the risetime of the signal . The noise bandwidth may then be expressed 
as 

K 1 

NBW = — — 
t ro 

Computer plots of the delay line output as expressed by Equation 6-11 

2 

show the signals to be approximately Cos in shape. If one assumes that 
a system bandwidth sufficiently wide to encompass 90% of the signal 
energy is required in order to faithfully reproduce it, then K-^ = 0.6. 

The output noise power in this bandwidth is therefore simply 

e M ■ G < e A *1* f* 2K ™° } 4 (Volts 2 ) 

The peak output signal amplitude is 

~ GR 0 ~ 

6 o = ~2~ X T * 

Defining the normalized peak output current for a given t^ as 


td U _ d T ( t dl) 
* l j J 


1 t (t dr°> 


we can write that (from Equation 6-3) 


? Q -1 tdK 
i T = r e f(— ) 

IT T 


The output peak-signal to rms noise voltage ratio is then 
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Using the simple model of Figure 6-8 for the effect of noise upon the 
timing accuracy of the CFD yields an estimate of the rms timing jitter: 

(1 - f.) t 
v 1 J ro 

0.8 S/N 


& 


where f^ if the fraction. 


( 2fl- fl )T tK Itro) 1/2 te A ,i A -s. .2KTR 0 ) 

0.8e~lQ tdi R 0 

f (— ) 0 


The time jitter described by Equation (6-14) is dependent upon the 
location of the peak of the charge distribution on the line. Except 
for the special case of an event centered on the line, it is necessary 
to evaluate the two delay line outputs and sum in quadrature the CFD time 
jitters to arrive at the total system jitter. For the purposes of simpli- 
city, we will assume that the event occurred at the center of the line and, 
hence, the overall jitter is given by Equation (6-14) multipled by vT. 
Computer studies of the delay line readout system indicate that the re- 
solution achieved with a charge located above the center of the chamber 
is maintained out to approximately a gap width from the edge of the active 
area where it begins to degrade due to the truncation of the mirrored 


cathode charge distribution. 
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FIGURE 6-8. 


Simplified model for calculating time jitter of constant 
fraction discriminator due to noise 
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The spatial resolution, a , achieved is given by 


1/2 


a = - — 
s tdi 


4 -) 


fl-fl)S 1/2 (tro/ T ) 

4.805 1 .■ (K lT ) A/ ^ ^ — 


1/2 


R 1/2 

T 7-^2° + 2KTR n ) 
(e A 2 + l A 4 0 


f t^i) 


R 


(6-15) 


for an event located at the center of the chamber. 


The spatial resolution clearly has a minimum (i.e. a minimum of the 
spatial inaccuracy) at the maximum of the function, 


g = 





1/2 


This function is plotted in Figure 6-9 for a 50 cm chamber with cathode 
wire spacing of 2 mm and an anode-cathode spacing of 5 mm. The model 
pulse time constant was 20 ns and the loading effect was included. 

The plot of Figure 6-9 indicates that there exists an optimum interwire 
delay of approximately 4 ns for the cathode readout line for this parti- 
cular MWPC. It should also be noted that this optimum is fairly broad, 
hence some flexibility in component selection is available to make use of 
standard values, for example. 

The last term of Equation (6-15) is dependent upon preamplifier and 
termination resistor noise levels. The resolution improves as , 
hence as large as delay line impedance as practical is desired. The 
spatial resolution of the chamber described above with a delay line of 
10000 characteristic impedance, a..dqlay of 4.5 ns/section (i.e. 2.25 ns/ 

-> - 1 ' 2 2 l / 2 / — 

mm) and a preamplifier with (e^ ) = 5 nV//Hz and (i ) =0.3 pA//Hz 

is given by a s = 4.22 x 10'^ mm-coul/Q for f^ = 0.5, = 0.6, and the 

"loading" effect considered. 

Considering only the electronic system, it appears that a resolution of 

-13 

100 microns can be achieved with an anode charge of 4.2 x 10 coul 
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FIGURE 6-9. Optimum delay per section for chamber with ^=250, N 0 =125 
s = 2mm, L=5mm and for a model-pulse time constant, 
t = 20 ns loading effect considered 




corresponding to approximately 2.7 x 10 electrons dependent, of course, 
on Q r 

6.5 Experimental Work 

An experiment was performed with our standard chamber in order to test 
the validity of the analysis in predicting the achievable spatial re- 
solution. The chamber was operated with a flowing gas mixture of 90% 
Argon-10% CO2 at one atmosphere. 

The delay line was designed to utilize standard component values and was 
found to have a delay of 8.8 ns/section (i.e. 4.4 ns/mm) and a char- 
acteristic impedance of 825 ohms. 

The test setup of Figure 6-6 was used with a micrometer mounted, colli- 
mated Fe^ source. The voltage preamplifier exhibited an estimated rms 
voltage noise of 5 nVv'iHz’ and current noise of 0.3 pA/Hz\ 

The experimental results are plotted in Figure 6-10 as a function of 
chamber high voltage and anode charge, which represents the relative gas 
gain as explained in Section 3. As expected, the resolution improves 
with increasing high voltage until it becomes limited by chamber pro- 
perties. These sources of resolution degradation are discussed in 
Section 5 . 

According to the above analysis, applied to this chamber configuration, 
the experimental MWPC readout system should exhibit a resolution given 
by 

-14 

a = 3.86 x 10 mm-coul/Q 

for the source centered on the chamber and with loading considered. The 
same system with a noiseless preamplifier would have a resolution given 
by 

-14 

a = 1.80 x 10 mm-coul/Q 

under the same conditions. Based on these results (considering only the 
readout system as limiting the resolution), it should be possible to 
achieve a 100 micron resolution for 2.4 x 10^ electrons using a practical 
preamplifier. 
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For an Fe^ X-ray source (5.9 KeV) and normal gases, approximately 220 

primary electrons are created requiring a chamber gas gain of approxi- 
4 

mately 1.8 x 10 . (When one speaks of gas gain in this context, it 
must be remembered that the readout system is only responsive to the 
cathode charge pulse for a time of the order of the total line delay 
time. The charge contained in the cathode signals in this period is 
some 40 to 50 percent of the total anode charge (for a ^ 1 ps long line) . 
The total delay of the experimental line was approximately 200 ns and the 
anode charge was measured with a charge sensitive preamplifier followed 
by single integration, double differentiation shaping of equal 1 ys time 
constants. The difference in these two "integration times" (i.e., that 
of the anode signal processing electronics and the delay line electrical 
length of 'v 200 ns) should yield an error of about 25% in the estimate 
of the anode charge. Because of this discrepancy, the actual resolution 
is correspondingly better than indicated by the graph of Figure 6-10. 



7. SPATIAL RESOLUTION MEASUREMENTS 


Resolution measurements were carried out with our standard chamber which 
was outfitted with a lumped element delay line which is described in de- 
tail in Section 6. The measurements were carried out with a col- 

limated Fe-55 source and a collimated B-source (Ru-106/Rh-106) . The 
experimental setup for measurements with the B-source is shown in Figure 
7-1. Each end of the delay line was viewed by a low noise voltage 
amplifier which feeds into a constant fraction discriminator (CFD) . The 
individual CFD outputs were used as start and stop pulses for a time to 
amplitude converter (TAC) whose output was accumulated in a multichannel 
analyzer. The analyzer then displayed the position distribution. The 
width of the distribution represents the spatial resolution. 

Fe-55 Source: The Fe-55 source was deposited on a 25y diameter wire and 

collimated by an adjustable slit. It could be moved across the thin 
chamber window in accurately measured steps by a micrometer mounted 
source holder. The width of the x-ray beam in the chamber gap was calcu- 
lated to be 

o = 45y 

B-Source: The B-source which was used to study the effect of inclined 

tracks, was similar to the Fe-55 source in construction. However, even 
with a well collimated source (o<50y in our case) it is difficult to 
achieve a well collimated electron beam in the sensitive region of the 
chamber because of multiple scattering of the electrons. To reduce the 
amount of scattering we mounted the source in a vacuum housing which was 
separated from the MWPC gap by a 1-mil Mylar window, as shown in Figure 
7-2. The source could be rotated to measure the spatial resolution as a 
function of the angle of inclination. The average r.m.s. width of the 
electron beam in the 1 cm gap at 0° inclination was calculated to be about 
150y for 3 MeV electrons. This increased to about 200y for 2 MeV elect- 
rons . 

For our measurements with the B-source we placed a plastic sintillator 
behind the chamber to select only high energy electrons above about 2 
MeV. The PHA was gated by the discriminator output. 
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FIGURE 7-2. Source arrangement for collimated B-source 


7.1 Results with Fe-55 Source 


Figure 6-10 from the previous section shows the results of resolution 
measurements with a collimated Fe-55 source in 90% Ar + 10% CO 2 as a 
function of the chamber voltage. The corresponding gas gain curve is 
curve b in Figure 4-4. Plotted is the observed f.w.h.m. of the position 
distribution and the spatial resolution which is defined as 

f.w.h.m 
a = 2.35 

Also shown is the resolution of the electronic readout system predicted 
by the delay-line analysis and the estimated contributions to the 
position uncertainty due to electron diffusion, lateral range of the 
primary ionization and source width. 

The somewhat surprising results are that the measured resolution does 
not significantly increase with increasing signal to noise ratio but 
rather flattens out to a value of 120y, and that the measured resolution 
does not approach the estimated theoretical limit. Since the noise of 
the system according to the delay line analysis is insignificant at the 
higher voltages, this may indicate that the spread of the discharge along 
the anode wire could be a significant limiting effect, which compensates 
the improvement in signal to noise ratio. 

Position distributions taken at four different positions at 2100 volts 
are shown in Figure 7-3. The distance between points was on the average 
3.3 mm. Figure 7-4 shows the position resolution at 2000 volts. The 
PHA calibration was 2 7. 5 y/ channel in both cases. 

7.2 Results with Electrons 

Measurements with the collimated electron source were carried out to in- 
vestigate the effect of inclined tracks on spatial resolution. The 
measurements were done with 90% Ar + 10% CO 2 at 2300 volts. The cor- 
responding beam width at 0° inclination was calculated to be a = 200y. 

The source width increased with increasing inclination due to the longer 
pathlength of the electrons in the MWPC gas: 

= 200 
cos^/2© 
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FIGURE 7-3 


Position distributions for collimated Fe-55 source at 
2100 volts in 90% Ar + 10% CC>2 . Distance between peaks: 
3 .3 mm 



FIGURE 7-4. Position distribution at 2000 volts 
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The source width was unfolded from the measured width of the position 
distribution according to 

2 2 2 

measured real source 

The measured and real resolution as a function of angle of inclination 
is shown in Table VII-1, which shows that the spatial resolution deter- 
iorates significantly with increasing angle of inclination. 

TABLE VII-1 

Resolution as a function of angle of inclination HV = 2300 volts, 90% Ar 

+ 10% co 2 . 


Inclination 

a meas . 

a real 

0° 

246y 

140y 

5° 

293y 

213y 

10° 

30Sy 

227y 

15° 

328y 

245y 

20° 

422y 

350y 

24° 

504y 

437y 


8. DELTA-RAY CALCULATIONS 
8 . 1 Introduction 

One of the limitations on the spacial resolution of a MWPC is the spatial 
distribution of the primary ionization produced by the passage of a 
charged particle through the MWPC gas. Specifically, secondary electrons 
(delta-rays) are produced which have energies such that their range is 
greater than the magnitude of the uncertainties in position location that 
we have been discussing. On the average, the delta-rays are produced 
symmetrically with respect to the particle trajectory and if the position 
determination was based upon the centroid of the charge collected in the 
MWPC, no systematic bias would be produced. However, the frequency of 
high-energy delta-ray production is low enough that the statistical 
fluctuations in delta-ray production might significantly increase the 
uncertainty of the position measurement of individual events. In addition. 
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if the MWPC does not linearly multiply the ionization deposited, due to 
space charge effects, the problem may be aggrevated. In effect, the ion 
pairs produced by delta-rays outside of the main core of primary ioni- 
zation would be more heavily weighted because of the reduction in multi- 
plication of the primary ionization column by space charge effects. The 
extreme case would be saturated chamber operation (e.g., "magic" gas 
operation) for which the MWPC response to an individual delta-ray would 
be equal to that for the primary ionization column. Finally, the whole 
problem is complicated if there is a magnetic field present. The Lorentz 
forces on the delta-rays are sufficient to radically alter their tra- 
jectories in a systematic way and to produce errors in the position deter 
mination. 

In order to make a quantitative estimate of the potential severity of 
these delta-ray effects, a computer program was generated that uses Monte 
Carlo techniques to perform a computer simulation of delta-ray effects. 
For this first-order assessment of the problem, a number of simplifi- 
cations were used to limit the complexity of the problem. The most sig- 
nificant of these simplifications are the neglect of the effects of 
multiple coulomb scattering and magnetic fields on the delta-ray tra- 
jectories. Delta-rays were propogated in a straight line from their 
point of production to the end of their range. In addition, range strag- 
gling was also neglected. Magnetic-field and multiple-Coulomb-scattering 
effects can be added to the computer simulation if necessary, but at the 
cost of a considerable increase in program complexity and computing time. 
The following is a description of the first-order computer simulation and 
a summary of the results for a simple detector configuration. 

8.2 Description of the Computer Program 

The geometry of the problem used in the computer simulation is shown in 
Figure 8-1. The formulae used in the computer simulation for delta-ray 
production by heavy, charged particles are given in Table VIII-1. The 
program- is set up to handle up to five layers of material along the inci- 
dent particle direction. The first material layer is treated as the 
sensitive region of the MWPC. The program simulates the passage of inci- 
dent particles through this combination of materials on an event -by-event 
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FIGURE 8-1. Geometry for delta-ray computer program 
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TABLE VIII-1. BASIC FORMULAE FOR DELTA-RAY PRODUCTION 


Definition of Symbols 

E : Incident particle kinetic energy per nucleon, 

z : Incident particle charge. 

8 : Incident particle velocity/velocity of light. 

Y : Incident particle total energy/rest mass. 

M : Atomic mass unit ( = 931.16 MeV) . 

m : Electron mass ( = .511 MeV). 

Z : Absorber atomic number. 

A : Absorber atomic weight. 

I : Absorber mean ionization potential, 

a) : Delta-ray kinetic energy. 

2 2 

W : Maximum delta-ray kinetic energy ( =2mB y .) 

oj o : Minimum delta-ray kinetic energy. 

0 : Angle of delta-ray emission with respect to incident particle 

direction. 

2 

R : Delta-ray empirical range (g/cm ) 


The probability that a heavy (mass »m) charged incident particle will 

produce a delta-ray within the energy interval of w to u + du while tra- 

2 

versing one gram/ cm of absorber is: 


d 2 n 

dxda> 


•15353 z 2 Z , , .2 

-Fa*- c 1 - 6 



The number of delta-rays with energy > co 0 


per g/cm 


2 . 
is : 


(8-1) 


dn i 
dx 


<»>- 


•15353, 2 Z r . 

T A 1 1 

8 2 A u o 



to 0 

w 


In 


O’ 


(8-2) 


2 

The ionization energy loss per g/cm due to collisions that produce delta- 
rays with energies < oj 0 is: 


dE 

dx 


) 


tor 


•15353z Z x , , Wco 0 > 0 2 x 

8 2 A ( 1" C - 6 > 


(8-3) 


The usual additional term for the density effect is neglected here because 
it is very small for practically all cases of interest. 
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TABLE VIII-1 (Continued) 


The relationship between the production angle and energy for delta-rays 
is 


cos 0 = - 5 - {u)/(oj + 2 m)} 
p 


1/2 


(8-4) 


The following empirical range-energy relationships are used for delta- 
rays (low energy electrons) . 


R = 


. 10 (1.265 - .0954 In u) 
.412a) 

.526(0 - .106 


( 0 < 3 MeV 

( 0 > 3 MeV 


(8-5) 
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basis. The spatial distribution of the energy deposited in the form of 
ionization in the MWPC sensitive region by delta-rays produced in all 
layers is calculated for each event and various distributions are ac- 
cumulated in the form of histograms for a specified number of incident 
particles (events) . 

The production of delta-rays by each incident particle is treated by 
starting with the incident particle at one boundary of the combination 
of materials. The delta-ray production coefficient (eqn. 8-2) is calcu- 
lated and the distance to the first collision which produces a delta-ray 
with energy above <o 0 is randomly generated from the corresponding ex- 
ponential distribution. If this distance is beyond the boundary for the 

layer in question, the process is reinitiated at the boundary of the next 

_2 

layer. Next, the delta-ray energy, m , is randomly generated from a to 
distribution which is a simplifying approximation to the actual energy 
distribution (eqn. 8-1). In turn, the angle of emission is calculated 
from eqn. 8-4, the range of the delta-ray is calculated from eqn. 8-5, 
and the projected range is checked to see if the delta-ray is able to 
penetrate to the first layer. If the delta-ray stops before reaching 
the first layer, it is neglected and the calculation starts again with 
generation of the point of the next delta-ray production. If the delta- 
ray can penetrate the first layer, the azimuth angle of emission is 
randomly generated from a uniform distribution and the trajectory is pro- 
jected into the x-y plane to determine the y-coordinate at entry into 
the first region. The range-energy relationships (eqn. 8-5) and the 
corresponding inverse energy-range relationships are used to determine 
the energy at entrance. The delta-ray is further propogated along a 
straight-line until it exists the first region or stops. The y-coordinate 
and energy at exit or the y-coordinate of the stopping point are calcu- 
lated. The energy deposition in the first region is the difference 
between the entrance and exit energies and the lateral position of the 
energy deposition is approximated by the y-coordinate of the mid-point 
of the delta-ray trajectory in the first layer. This sequence is repeated 
until the incident particle has passed through the entire combination of 
materials. The several possibilities depending upon whether the delta-ray 
was produced in the first layer or prior to the first layer and whether 
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the delta-ray stops in the first layer or exists the first layer are 
graphically depicted in Figure 8-1. All lateral positions are projected 
into one plane because this is the quantity measured by the MWPC systems 
with which we are dealing. 

The ionization energy loss in the first layer due to incident particle 
collisions in which the energy transferred to the secondary electron is 
less than to 0 (so called soft collisions) is calculated according to eqn. 
8-3. This energy deposit is treated as if it is concentrated along the 
incident particle trajectory (i.e., y = 0) . The MWPC is assumed to re- 
spond linearly to the ionization produced by delta-rays and the MWPC gain 
nonlinearity is approximated by assuming that the MWPC response to the 
soft-collision ionization is proportional to the square-root of the 
energy deposited, in correspondence with our experimental results (see 
Section 9) . 

For each incident particle, the lateral centroid of the energy deposition 
in the first layer by delta-rays produced in each of the layers and for 
the energy deposition from all layers is calculated. This distribution 
of lateral centroids of energy deposition is the primary result of the 
computer simulation because it represents the distribution of position 
measurements which a true centroid-sensitive MWPC system would produce. 

In addition, the lateral distributions of energy deposition in the first 
layer due only to delta-rays produced in each of the layers, and the 
distributions of the fraction of the total energy loss in the first layer 
due to delta rays produced in each of the layers are also generated. 

These distributions are mainly useful for diagnostic purposes and for 
checking the relative importance of the individual layers. 

8.3 Results of the Computer Simulation 

The computer program described above was run for the following simple 
configuration: a 1-cm. thick layer of Argon (the MWPC sensitive region) 
preceded by a 0. 001-in thick Mylar layer (the MWPC window) and a 1-m. 
thick layer of air. One-thousand event runs for incident particles with 
E = 2.1 GeV/nucleon and z = 4, 8, 12, 16, 20 and 26 were performed. A 
summary of the results is contained in Table VIII-2 and representative 
distributions generated are shown in Figures 8-2 through 8-5. 
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ENERGY LOSS/(INCIDENT PARTICLE CHARGE) (KeV/mm) 
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FIGURE 8-3 

PROJECTED LATERAL DISTRIBUTION OF THE CENTROID 
OF THE ENERGY DEPOSITION IN A 1 - cm THICK 
ARGON LAYER FOR Z = 4, 8, 16, AND 26. THE 
CENTROID IS CALCULATED WITH SQUARE-ROOT MWPC 
RESPONSE TO THE SOFT COLLISION IONIZATION AND 
A LINEAR RESPONSE TO DELTA-RAYS PRODUCED IN 
THE ARGON. 


INCIDENT ENERGY = 2.1 GeV/NUCLEON 
MINIMUM DELTA-RAY ENERGY = 10 KeV 
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ENERGY LOSS/(INCID£NT PARTICLE CHARGE ) 2 KeV/c 
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NUMBER OF EVENTS 
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PROJECTED LATERAL DISTRIBUTION OF THE CENTROID OF 
THE ENERGY DEPOSITION IN A 1 - cm THICK ARGON 
LAYER FOR Z = 4, 8, 16, AND 26. THE CENTROID IS 
CALCULATED WITH SQUARE - ROOT MWPC RESPONSE TO 
THE SOFT COLLISION IONIZATION AND A LINEAR 
RESPONSE TO DELTA-RAYS PRODUCED IN THE ARGON 
AND IN A 1 - mil THICK MYLAR WINDOW AND 1 -m 
THICK AIR LAYER PRECEDING THE ARGON. 

INCIDENT ENERGY = 2.1 GeV/NUCLEON 
MINIMUM DELTA-RAY ENERGY = 10 KeV 
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The figures showing the projected lateral distributions of the energy 
deposition in the MWPC sensitive region (layer 1) due to delta-rays pro- 
duced in layer 1 and layer 3 (the 1-m air layer preceding the MWPC) are 
valid for all incident charges. For any value of incident charge, z, 
the ordinate should be multiplied by z . No distribution is given for 
layer 2 (the Mylar window) because its contribution was negligible in all 
cases. In order to get a feel for the contribution to the total energy 
deposition due to delta-rays relative to the energy deposition due to soft 
collisions, it should be noted that the energy deposition due to soft col- 
lisions (energy transfers < a> 0 = 10 KeV) divided by z is equal to 1.87 
KeV. Delta-rays with an energy < 10 KeV are emitted at angles > 80 de- 
grees and have a range < 1mm in Argon. Therefore, the energy deposition 
due to soft collisions is contained in a column of radius < 1mm and has 
a mean projected lateral energy deposition density divided by z 2 equal 
to .935 KeV/mm, or about 20 times the maximum density due to delta-rays. 
Although the lateral energy deposition density due to delta-rays is much 
lower than that due to soft collisions, the delta-ray contribution is 
spread over very large lateral distances. This is especially true of the 
contribution from the 1-m air layer as would be expected from the geometry 
of the problem. When we consider the lateral distribution of the centroid 
of the energy deposition, we see that the delta-rays from the air layer 
produce a long tail on the lateral distributions. Our treatment of the 
non linearity of the MWPC response produces an accentuation of the delta- 
ray effects because while the delta-ray contributions are increasing as 
z 2 as z increases, the soft collision contribution is only increasing as 
z. However, even if the MWPC responded linearly the delta-rays will pro- 
duce a position measurement distribution which is definitely not Gaussian. 
This can be seen by comparing the rms and 68% widths for low z in Table 
VIII-2. If the distributions were Gaussian the two widths would be equal. 
The fact that 68% widths are smaller than rms widths reflects the fact that 
the distributions have a central spike with very broad, low-level tails. 

The fact that the widths of the centroid distributions level -off or even 
decrease as z increases above 12 to 16 comes about because the delta-ray 
production is becoming frequent enough that statistical fluctuations are 
starting to average out in the centroid. 
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TABLE VI I I -2 


Configuration Definition 


Layer 

Material 

Density 

(g/cc) 

Thickness 

(cm) 

1 

Argon 

1.784 x 10' 3 

1.0 

2 

Mylar 

1.18 

2.54 x 10" 3 

3 

Air 

1.205 x 10- 3 

100.0 


Incident Particle Definition 

Energy = 2100 MeV/nucleon 
Maximum delta-ray energy = 9.81 MeV 
Minimum delta-ray energy = 0.01 MeV 


Charge 



Units 

4 

8 

12 

16 

20 

26 

Mean soft collision energy de- 
posited in layer 1 

KeV 

29.9 

119.6 

269.1 

478.4 

748 

1263 

Mean energy deposited in layer 
1 by delta-rays produced in: 

KeV 







layer 1 


4.2 

15.3 

36.9 

59.1 

95.1 

161.8 

layer 3 


1.9 

7.1 

16.3 

29.2 

45.0 

78.0 

Mean total energy deposited in 
layer 1 

KeV 

36.0 

142.0 

322.3 

566.7 

888 

1503 

Widths of projected lateral 
distribution of energy depo- 
sition in layer 1 by delta- 
rays produced in: 








layer 1 - rms 

cm 

.44 

.45 

.43 

.44 

.43 

.43 

- 68%* 


.34 

.28 

.32 

.31 

.28 

.30 

layer 3 - rms 


41.7 

42.9 

42.0 

42.8 

43.1 

43.5 

- 68% * 


40 

42 

42 

42 

42 

42 

Widths of projected lateral 
distribution of the centroid 
of energy deposition in layer 
1 including delta-rays pro- 
duced in: 








layer 1 - rms 

cm 

.12 

.15 

.15 

.15 

.13 

.11 

only - 68%* 


<.02 

.04 

.08 

.09 

.10 

.10 

all layers - rms 


8.63 

9.69 

7.17 

6.66 

5.38 

4.53 

- 68% * 


<.2 

5.0 

5.3 

5.4 

4.4 

3.9 


+y 

* 68% width = y such that / contains 68% of the events. 

-y 
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As previously noted, the effects due to delta-rays produced in materials 
such as the air layer preceding the MWPC may be greatly reduced by the 
presence of a magnetic field component perpendicular to the incident 
direction which prevents delta-rays from moving very far in the incident 
particle direction. Also many of the delta-rays produced in the air 
layer penetrate the MWPC at distances which are far enough from the point 
at which the incident particle traverses the MWPC to produce separated 
pulses on the delay line. In this case of separated pulses, it is pos- 
sible that the effects of the delta rays could be electronically minimized. 
However, even if we limit ourselves to the effects of the delta-rays pro- 
duced in the MWPC gas itself, the delta-rays produce a rms position un- 
certainty of about 1 to 1.5 mm for a 1 cm thick sensitive region. If 
we consider the 68 % width as a better indication of the spacial resolution, 
we see that the delta-ray effects are essentially the order of or less 
than other limits discussed in the section on the spacial resolution for 
z 4, and rise to about 1 mm for z = 26. 

8.4 Energy Loss Calculations for Gas Mixtures 

In order to calculate dE/dX for the mixture we need the mean ionization 
potential and Z/A of the mixture. If the fraction by weight of each ele- 
ment in the mixture is denoted by f^ 

Z , _ _ 

A ^mix £ i A j 


In = I f i C|) In Ii/Z f 


mix 


With I . , the dE/dX is determined by looking up the stopping power for 
mix 24 

*mix a P roton with equal velocity in Barkas § Berger , Table I and 

multiplying by Z. } 

A mix. 


For nuclei with 2.1 GeV/nucleon kinetic energy the equivalent (i.e. equal 
velocity) proton energy is 2.12 GeV. 


We have worked 

with the 

following mixtures: 




Constituent 

Ar 2 

C0 2 C 4 H io 

X e 




Mix 






1 

0.94 

0.06 

- 




2 

0.90 

0.10 

- 




3 

0.80 

0.20 

- 

(fraction by volume) 

4 

0.90 

0.10 

- 




5 

0.75 

0.25 

- 




6 

0.81 

0.05 

0.14 




The four gases 

in the mixture have the following properties 




Ar 

co 2 

c 4 h 10 


Xe 

Density 

(mg/cc) 

1.784 

1.977 

2.673 

4 

.587 

Molecular Wgt. 

Cg) 

39.95 

44.01 

58.12 

131.30 

The mixtures have the following properties: 




Constituent 

Ar 

co 2 

C 4 H 10 

Xe 

Density of 

Mix 





Mix (mg/cc) 

1 Partial 

Density 

1.677 

.119 

- 

- 

1.796 

(Fraction by Weight) (.934) 

(.066) 




2 


1.606 

.198 

- 

- 

1.804 



(.890) 

(.110) 




3 


1.428 

.385 

- 

- 

1.813 



(.786) 

(.214) 




4 


1.606 

- 

.267 

- 

1.873 



(.857) 


(.143) 



5 


1.339 

- 

.668 

- 

2.007 



(.665) 


(.335) 



6 


1.440 

.103 

- 

.643 

2.186 



(.659) 

(.047) 

( 

.294) 
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Finally we have for the mixtures in terms of fractional weights by ele- 
ments : 



Element H 

C 

0 

Ar 

Xe 

Mix 






1 

- 

.018 

.048 

.934 

- 

2 

- 

.030 

.080 

.890 

- 

3 

- 

.058 

.156 

.786 

- 

4 

.025 

.118 

- 

.857 

- 

5 

.058 

.277 

- 

.665 

- 

6 

- 

.013 

.034 

.659 

.294 

Z/A 

1.0 

.5 

.5 

.451 

.411 

I 

(eV) 19 

78 

103 

210 

555 

The 

results of the 

equations 

above are: 



Mix 

Z/A . 
mix 

I . 
mix 

Stopping 

power 

dE MeV 

dX ^ g/cm^ 


1 

.454 

199 

3.357 

1.521 


2 

.457 

188 

3.373 

1.540 


3 

.462 

175 

3.400 

1.568 


4 

.471 

161 

3.432 

1.614 


5 

.497 

119 

3.530 

1.752 


6 

.442 

288 

3.235 

1.430 



These are mean loss rates. A correction needs to be applied for escaping 
delta-rays and entering delta rays from materials preceding the MWPC. 

For pure Argon = 1.505 MeV/g/cm^ 

and hence AE in 1 cm = 2.68 KeV 

2 

The delta-ray calculations gave an energy loss /Z =2.23 KeV 
Therefore our best estimate of energy deposited is .832 of mean dE/dX. 

To calculate the no. of ion pairs produced in the MWPC gas we also need 

to know the mean energy required to produce an ion pair. The values in 

25 

the summary table below were calculated following Hurst et al_ : 
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Gas Mixture 


94% Ar 90% Ar 80% Ar 90% Ar 75% Ar 81% Ar 

6% C0 2 10% C0 2 20% C0 2 10% 25% 14% Xe 






C 4 H 10 

C 4 H 10 

5% C0 ; 

Density (mg/cc) 

1.796 

1.804 

1.813 

1.873 

2.007 

2.186 

Mean dE/dX (MeV/ 
g/cm^) for Z=l, 

B = .952 

1.521 

1.540 

1.568 

1.614 

1.752 

1.430 

Energy Deposi- 1.265 
tion Rate (MeV/ 
g/cm 2 ) dE/dX .0.832 

1.280 

1.302 

1.341 

1.458 

1.189 

Energy Deposited 
per cm @ STP 
(KeV/cm) 

2.27 

2.31 

2.36 

2.52 

2.93 

2.60 

Mean energy to 
form an ion 
pair (eV) 

26.2 

26.4 

27.2 

24.3 

25.4 

24.0 

Ion pairs per 
cm (cm“l) for 

87 

88 

87 

104 

115 

108 


Z=l, 6 = .952 

9. BEVATRON MEASUREMENTS AND RESULTS 

A series of measurements was carried out with several of our standard 
MWPC's in a relativistic heavy ion beam at the Bevatron. Two heavy ion 
beams were available at the Bevatron: 

2.1 GeV/nucleon C ^ 
and 

2.1 GeV/nucleon 0 16 

The purpose of our measurements was to measure the gas gain curves for 
12 16 

C and 0 nuclei in various gases for comparison with the gas gain 
curves obtained with Fe-55 x-rays. We also wanted to study the pulse shapes 
of the MWPC current pulses, charge pulses and delay-line pulses as a 
function of the MWPC gas and supply voltage. 

The experimental setup contained 4 chambers as shown in Figure 9-1. The 
chambers were mounted on a common base plate with electrical and gas feed 
throughs. The base plate itself was mounted to a turn-table. A hood made 
out of aluminum foil was placed over the setup to act as an RF-shield. 
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FIGURE 9-1 


Experimental setup for Bevatron measurements 
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Three of the chambers (chambers No. 2, 3, and 4 in Figure 9-1) were out- 
fitted with a capacitively coupled, bulk delay line of the Perez-Mendez 
type and chamber No. 1 was connected to a directly coupled, lumped ele- 
ment delay-line described in Section 6. The anodes of chamber 2, 3 and 
4 were connected to Ortec 109-PC charge sensitive amplifiers. The cathode 
of each of the chambers 2, 3 and 4 which was not coupled to the bulk de- 
lay-line was connected to a fast current amplifier. A current amplifier 
was also connected to the anode wires of chamber No. 1. Figure 9-2 shows 
the block diagram of the experimental setup. A coincidence telescope, 
which consisted of two 1 cm x 1 cm plastic scintillators, defined the 
particle beam and was used to trigger the oscilloscope camera and to gate 
the pulse height analyzer. 

The gas gain was measured in three different gases, 

80% Ar + 20% C0 2 

75% Ar + 25% Isobutane (C^^g) 

81% Ar + 14% Xe + 5% C0 2 (for C 12 - runs 
only) 

with C 12 and 0 ^ nuclei and Fe-55 x-rays. 

The pulse shapes of the current pulses, delay-line pulses and charge 
pulses as a function of HV were photographed for C* 2 and 016 nuclei in 
80% Ar + 20% C0 2 and 75% Ar + 25% Isobutane, and for C* 2 nuclei only in 
81% Ar + 14% Xe + 5% C0 2 . 

The results of the gas gain measurements are plotted in Figure 9-3 for 
C 12 and in Figure 9-4 for O 1 ^. The ordinate of these two figures is the 
number of electrons collected at the MWPC anode as measured with the 
amplifier system using 1 ysec shaping time constants. These data can be 
converted to relative* gas gain by dividing the number of electrons col- 
lected by the number of primary electron-ion pairs produced in the MWPC 
sensitive region (see Section 8.4). In Figure 9-5, we show the relative 
gas gain for Fe-55, C* 2 , and 0^ in 80% Argon and 20% C0 2 . We see, as 
can be expected from our considerations of space charge effects in Section 
4.4, that the onset of deviation from strict proportionality occurs at 
lower gas gains with increasing primary ionization densities. The 

*Recall that, as discussed in Section 3, an additional correction for 

the fraction of the total charge collected in the ' 1 ysec time constant 

of the amplifier would have to be made to obtain absolute gas gains. 

96 



DELAY-LINE 



FIGURE 9-2. Block diagram of Bevatron experimental setup 
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12 

FIGURE 9-3. Gas gain for 2.1 GeV/nucleon C -nuclei 
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Figure 9-4. 


Gas gain for 2.1 GeV/nucleon 0 


16 


nuclei 




RELATIVE GAS GAIN 



FIGURE 9-5. 


12 1 A 

Comparison of gas gain of Fe-55 x-rays and C and 0 
nuclei in 80% Ar + 20% CO 2 
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magnitude of the reduction in gain was proportional to the total charge 
in the signal. At a given high voltage, the charge in the signal is 
proportional to the primary ionization, and hence the reduction in gain 
is an increasing function of the primary ionization density. We also 
see that the relative reduction in gain increases as the voltage increases, 
again as expected, since the amount of charge in the signal increases with 
increasing voltage. 

Another way of looking at the saturation characteristics of our chambers 
is to plot the measured number of electrons for Fe-55, C ^ an d 0^ for 
the three gases as a function of the energy loss in the chamber gap. This 
was done in Figure 9-6 on log-log paper. The energy loss is expressed in 
terms of the energy loss of a minimum ionizing particle and is corrected 
for delta-ray effects, (see Section 8.4). Strict proportionality would 
be represented by a slope of + 1 in this plot. The effect of saturation 
is to reduce the slope of the InN vs InAE plot. In Figure 9-7 we have 
plotted the slopes from Figure 9-6 as a function of the number of electrons 
collected for a minimum ionizing (z = 1) particle. The increasing 
saturation with increasing gain shows very clearly and is seen to be ap- 
proximately equal for the various gases. We have already noted that 
this deviation from proportionality will accentuate the detrimental ef- 
fects of delta-rays on spatial resolution. It is, therefore, desirable 
to operate MWPC's at the lowest possible gain that still affords enough 
signal to achieve the required spatial resolution for the least ionizing 
particles of interest. We have seen in Section 6.5 that a signal of ap- 
proximately 2 x 10^ electrons is required for a spatial resolution ap- 
proaching 100 pm with the lumped element delay line. If we require that 
signal for minimum ionizing particles, we see from Figure 9-7 that the 

MWPC response as a function of primary energy deposition in the MWPC will 

1/2 

be approximately AE (i.e. slope - 0.5). This square-root MWPC response 
was what we used in the analysis of the delta-ray problem in Section 8. 


NUMBER OF ELECTRONS 



AE (X MIN. IONIZING) 


FIGURE 9-6. Measured number of electrons v.s. energy loss of Fe-55 x- 
rays and C 12 - and 0 16 -nuclei in three different gases as a 
function of HV. 
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In Section 2 we noted that the desired range of incident charge for 
cosmic ray spatial detectors is commonly 1 <_ z <_26. Our Bevatron measure- 
ments obviously do not answer the question of whether the response ob- 
served for 1 < z < 8 can be extrapolated to z = 26. Recall that in 
Section 4.5 we found that the space charge saturation limit on the signal 

9 

charge density was about 10 electrons/cm (i.e., about 0.5 esu/cm) for 
our MWPC electrode geometry with operating voltages of about 2200V. Since 
we expect the effective width of the discharge along the anode to be the 

O 

order of 0.1 cm, we should reach space charge saturation at about 10 
electrons. We can see that all of the curves in Figure 9-6 extrapolate 

g 

to less than 10 electrons for a AE value corresponding to z = 26, and 
hence such an extrapolation at least does not violate the space charge 
saturation limit. 

10. MWPC - DESIGN CONSIDERATIONS 


In Section 2 we discussed some of the requirements for a typical spatial 
detector system for a high energy cosmic ray experiment, which in this 
case was the superconducting magnetic spectrometer experiment on HEAO-B. 

One of the important requirements for the MWPC frame in that experiment 
was the minimum septum width of the frame along the wall of the cryostat. 
The frame material had to be non-magnetic and had to exhibit low outgasing. 


The minimum septum width for a given frame material is dependent on the 
frame load due to the tensioning of the anode and cathode wires. To calcu- 
late the frame load due to the anode wire tensioning one has to consider 
the electrostatic problems in MWPC's, i.e. the electrostatic stability 
of the anode wires. The minimum tension, t c , required to keep the anode 
wires in equilibrium is 


t 


c 



( 10 - 1 ) 


where q is the charge per unit length, 1 is the wire length and s the 
anode wire spacing. 

The charge per unit length was given previously as 


q = 


Ya 

2X 
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V q is the chamber voltage and X is a function of the geometrical configu- 
ration: 


X = 



In 


trd 

s 


where L is the anode to cathode spacing and d is the anode wire diameter 


The total frame load with wires at the critical tension is 




( 10 - 2 ) 


To calculate the frame load for a specific case we choose the large area 

chamber of Section 2 with 1 = 66 cm, with our geometrical configuration, 

-3 -3 

i.e. L = 0.5 cm, s = 0.4 cm and d = 20 x 10 cm, 50 x 10 cm and 100 x 
IQ" 3 cm. 


The critical tension and the total frame load are plotted in Figure 10-1 
as a function of the operating voltage V Q which determines the charge 
per unit length. The solid curves cover the voltage range for which the 
relative gas gain was measured for Fe-55 x-rays in 80% Ar + 20% C0 2 . 

(See Figure 4-2) . The arrow indicates at which point the relative gas 

3 

gain was 3.5 x 10 . From the results of Figure 10-1 and equation (10-2), 
it is obvious that it is desirable to use small diameter, widely spaced 
anode wires to minimize the total frame load. According to our earlier 
results this condition is also preferable from an operational point of 
view, i.e. small diameter, widely spaced anode wires result in an in- 
crease in the maximum charge per unit length that can be accumulated and 
thus increases the strictly proportional region of the chamber. 

The actual tension that will be put on the anode wires is above the 
critical tension but below the yield tension, so that any changes in 
wire tension due to temperature effects will keep the tension in a safe 
range from an operational and mechanical point of view. 

Tensioning of the cathode wires is normally not as critical as tensioning 
of the anode wires. The tension mainly has to be large enough to prevent 
sagging of the cathode wires. In MWPC's with integral delay-line readout, 
however, tensioning of the cathode-or sense wires can be a problem since 
wires of rather large diameter (e.g. 150y) have to be used to minimize 
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FIGURE 10-1 Critical anode wire tension t c (in grams) and total frame 
load (in pounds) for a large area MWPC as a function of 
HV, arrows indicate HV at which the relative gas gain for 
Fe-55 x-rays is 3.5 x 10 3 







the attenuation of the delay-line. Considerable tension has to be put on 
large diameter wires simply to keep them straight, and this can lead to 
rather large frameloads. 

The frame material has to be limited to nonmagnetic materials, and since 
the frame also has to insulate anode from cathode planes, G-10 fiberglass 
is usually selected for the frame material. The disadvantages of G-10 are 
its outgassing characteristics and the possibility of creeping under long 
term, high loading. Other materials which are suitable are Mycalex-400, 
which is used as frame material for large area spark chambers for satellite 
experiments, and glass. Glass has very good mechanical and thermal prop- 
erties and is very low in outgassing. Table X-l summarizes some of the 
properties of G-10, Mycalex-400, and glass. 
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